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Institute for Thermal Energy Technology and Safety 


Structure and Activities of the Institute for Thermal Energy 


Technology and Safety (ITES) 


Walter Tromm 


The Institute for Thermal Energy and Safety 
(ITES) is situated with its offices and research 
laboratories in Campus North of Karlsruhe In- 
stitute of Technology (KIT). 


The research of ITES focuses on conversion 
from thermal power to electric power for future 
power systems without CO2 emission includ- 
ing hydrogen technologies and on safety. In 
the past ITES concentrated on safety features 
and on methods to mitigate severe accidents 
for nuclear power plants. Nowadays, ITES 
uses these competences, gained from a pro- 
found experience in numerical simulation and 
in design and operation of large-scale experi- 
mental facilities to apply this knowledge mainly 
in the area of renewable energies. 


The combination of science and technology 
with education and training is a systematic ap- 
proach at KIT, and ITES is contributing accord- 
ingly to courses in mechanical engineering, su- 
pervises several bachelor and master theses 
each year and coordinates master programs in 
energy technologies. Compact courses on en- 
ergy technologies are given also in executive 
master programs and in the Framatome Pro- 
fessional School, which is funded by industry 
and managed by ITES. 


Resources 


The employees of ITES contribute to the re- 
search-programs Renewable Energies (EE), 


Employees assigned according to 
Helmholtz programs, 2020 


gm Fusion 
m NUSAFE 
m EE 

m SCI 


Fig. 1 Assignment of ITES personnel to helmholtz-programms 
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Storage and  Cross-linked Infrastructures 
(SCI), Fusion (FUSION), and Nuclear Waste 
Management, Safety and Radiation Research 
(NUSAFE) in the research field energy of the 
Helmholtz Association (HGF). By the end of 
2020, around 80 scientists, engineers and 
technicians have been working at ITES on this 
wide range of CO»-free technologies for energy 
conversion. Fig. 1 illustrates that still about 60 
96 of the ITES employees were working for the 
NUSAFE program, but steadily reduced in re- 
cent years. Around two third of the employees 
were funded in 2020 by the Helmholtz Associ- 
ation (HGF), the others by third party funds of 
the European Commission, by industry, by 
German ministries or by other research funds. 
Doctoral students as well as students of the 
Baden-Wuerttemberg Cooperative State Uni- 
versity (DHBW) were filling around 2096 of 
these positions at ITES. In addition, students 
perform their bachelor or master theses or 
spend an internship in the research laborato- 
ries of ITES. 


An overview of the structure of ITES is given 
by the organization chart, Fig. 2. Because all 
working groups are embedded in the research 
field energy of the Helmholtz association 
(HGF) a close collaboration within the groups 
is guaranteed. 


Working groups on Severe Accident Research, 
on Multi-Phase Flow as well as the Framatome 
Professional School still concentrated in 2020 
primarily on nuclear applications, whereas the 
Karlsruhe Liquid Metal Laboratory (KALLA) 
and the Hydrogen group were addressing nu- 
clear as well as renewable energy technolo- 
gies. The working group on Magneto-Hydrody- 
namics is primarily working on nuclear fusion 
applications, whereas the working group on 
Energy and Process Engineering is concen- 
trating on geothermal energies. The working 
group Multiphase Flow is mainly experimen- 
tally investigating technically relevant, complex 
flow phenomena. The field of activity currently 
focuses on the build-up of the complex struc- 
ture of the COSMOS-H high-pressure loop. 
The working group on Accident Management 
Systems continued to extend the application of 


its simulation models to several critical infra- 
structures. Thus, the institute covers a wide 
field of different energy technologies and re- 
lated safety investigations, and the share of 
personnel resources on the different research 
topics reflects the requirements of the Helm- 
holtz Association's energy research field. 


Institute for Thermal Energy 
Technology and Safety 


W. Tromm (acting) 
Deputy W. Raskob 


Multiphase 
Flow 


S. Gabriel 


Severe Acci- 
dent Research 


X.-Y. Gaus-Liu 


Accident 
Management 
Systems 


W. Raskob 


Liquid Metal 
Laboratory 


T. Wetzel 


Framatome 
Professional 
School 


A. Class 


Magnetohy- 
drodynamics 


L. Bühler 


Hydrogen 
T. Jordan 


Energy and 
Process 
Engineering 

D. Kuhn 


Infrastructure 
K. Arheidt 


Fig. 2: Organization chart of the Institute for Thermal En- 
ergy Technology and Safety 
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They are all supported by a joint infrastructure, 
comprising a metal workshop, manufacturing 
urgent test components, a welding shop, and 
an electromechanical workshop. Other tasks of 
the infrastructure include the IT-administration, 
business administration and public website of 
the ITES. The Infrastructure team is active as 
well in education and training activities. 


Every year, at least 4 students of the Baden- 
Württemberg Cooperative State University are 
employed by ITES, managed by the Infrastruc- 
ture group, to work with the research teams as 
part of their educational program. 


Working in a research area with industrial ap- 
plication ITES practices a Quality Management 
System appropriate to EN ISO 9001. 


Group: Magnetohydrodynamic 


Magnetohydrodynamics for liquid-metal blankets 


L. Bühler, H.-J. Brinkmann, C. Courtessole, V. Klúber, C. Kóhly, C. Mistrangelo 


Introduction 


The design of breeding blankets represents a 
major challenge for fusion reactor engineering 
due to performance requirements and severe 
operating conditions in terms of heat load and 
neutron flux. Liquid metal alloys such as lead 
lithium PbLi are considered as breeder mate- 
rial due to their lithium content, as neutron mul- 
tiplier, and neutron shield. Because of their 
large thermal conductivity and the possibility to 
be operated at high temperature, liquid metals 
may serve additionally as coolants to remove 
the nuclear heat released in fusion reactors. 
On the other hand, when the electrically con- 
ducting fluid moves in the strong magnetic field 
that confines the plasma, electric currents are 
induced. The latter are responsible for strong 
electromagnetic Lorentz forces and high pres- 
sure drop. Magnetohydrodynamic (MHD) pres- 
sure drops have to be carefully quantified, 
since excessive values can jeopardize the fea- 
sibility of the considered blanket concept. 
Moreover, Lorentz forces affect the velocity 
distribution and suppress turbulence with impli- 
cations on heat and mass transfer. 


The physical parameters characterizing the 
MHD flow are the Hartmann number, the 
Reynolds number and the Grashof number 


L ATL 
Ha=BL p. docu qo AME 
pv V V 


where B, L and uo denote magnetic field, typical 
length scale and characteristic velocity. The 
fluid properties are density p , kinematic vis- 
cosity v , electric conductivity o and volumet- 
ric thermal expansion coefficient £. 


Ha? and Re stand for ratios of electromagnetic 
to viscous forces and inertia to viscous forces, 


respectively. The importance of buoyancy 
forces is quantified by Gr, where AT is a typi- 
cal temperature difference and g stands for 
gravity. 


The MHD work at ITES KIT supports DEMO 
blanket engineering activities as part of the 
EUROfusion consortium and the test blanket 
program of Fusion For Energy (FAE) for ITER. 


Numerical simulations for liquid metal 
MHD flows 


Simulation of 3D MHD flows in pipes with flow 
channel inserts 


When duct walls are thick and electrically con- 
ducting, flow-induced currents find shortcuts 
along these walls, resulting in high current den- 
sity in the fluid, strong electromagnetic forces 
and considerable pressure drop. In order to 
mitigate MHD effects in fusion blankets, elec- 
trically insulating flow channel inserts (FCI) 
have been proposed. They are loosely fitted in 
liquid metal ducts or pipes with the purpose of 
electrically decoupling the fluid region from the 
well-conducting blanket structure. [1]. In the 
following, a FCI of sandwich-type is consid- 
ered, as schematically shown in Figure 1, 
which consists of an insulation layer (orange), 
enclosed by two thin sheets of steel to protect 
the insulation from fluid infiltration. Of particular 
interest are 3D MHD effects present when the 
flow enters or exits the FCI or at gaps between 
two FCls. At these locations current density is 
locally increased by extra 3D currents j3D, 
which affect the flow. This leads to local modi- 
fications of electric potential, pressure and ve- 
locity. Numerical investigations of complex FCI 
entry flows at fusion relevant parameters 
Ha —2000, Re=20000, allows to quantify 
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these 3D effects and provide insight into re- 
lated experiments performed in the MEKKA la- 
boratory at KIT. 


Figure 1 Sketch of the MHD flow entering a FCI, as 
tested for a circular pipe in the MEKKA laboratory. The 
orange surface indicates the insulation layer. 


The computational model considers the FCI 
entrance flow at a high degree of geometrical 
detail, which allows the resolution of wall and 
conducting sheets of the FCI for determination 
of currents. The geometry, which takes into ac- 
count the liquid metal filled gap between FCI 
and pipe wall, is shown in Figure 1. 


In Figure 2 the electric potential distribution is 
shown in the lower half of the pipe. The incom- 
ing flow is fully developed far upstream with no 
variations of potential along x and along mag- 
netic field lines. Here, currents close in cross 
sectional planes, as shown by current lines dis- 
played in the upper half of the geometry at x=xı, 
x». Potential and currents in the other halves 


of the pipe behave in a y- symmetric way and 
they are in agreement with known solutions for 
fully developed MHD flow under a uniform 
magnetic field (e.g. [2]). 


After entering the FCI for x>0, the insulation 
prevents fluid currents from entering the thick- 
wall domain. They are blocked by the insula- 
tion and forced to close along the thin inner 
sheet of the FCI (red current lines). Far enough 
downstream, the flow approaches again fully 
developed conditions. 


Three dimensional effects near the FCI en- 
trance at x-0 lead to axial components of cur- 
rent density in the fluid and in the wall. This is 
caused on one hand by the fact that the poten- 
tial in the outer wall drops to zero for x>>1, 
since no currents are supplied to this part of 
the pipe due to the presence of the insulating 
layer. The resulting axial potential gradient 
drives 3D currents in the outer wall, as shown 
by the strongly bent blue current lines near x=0. 
On the other hand, current loops closing along 
the thin conducting sheet of the FCI experi- 
ence higher resistance and therefore the po- 
tential difference in the cross section is higher 
than in the well conducting part of the pipe for 
x«0. The reason is the significant difference of 
wall conductivity. 3D effects on currents and 
potential persist downstream into the FCI over 
multiple characteristic lengths. 


insulation 


Figure 2 MHD flow entering a FCI for Ha=2000, Re=20000. Electric potential distribution in the lower half of the pipe 
and electric current streamlines in the upper half. Current streamlines closing over the wall domain are displayed in 
blue, while those closing along the FCI are drawn in red. The orange surface represents the insulation layer and 


the black straight line at x=0 marks its beginning. 
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eddies 


Figure 3 3D profiles of velocity at three axial positions near the entrance of the FCI. The insulation starts at x=0. 


Such complex currents interact with the mag- 
netic field, resulting in Lorentz forces that affect 
the flow. Figure 3 shows the resulting velocity 
distribution across the FCI entry region at three 
axial positions downstream in the FCI. The in- 
coming flow, not shown in the figure, has a uni- 
form core velocity and thin boundary layers 
along the pipe wall. 


As soon as the fluid enters the insulated part 
(x0), it moves towards the sides forming a 
concave velocity profile along the transverse z- 
direction. This rapidly changes into a convex 
profile with a center peak (x-1). Then, over 
several characteristic lengths, the velocity 
slowly converges towards a fully developed 
flow in the FCI, where the core velocity is al- 
most constant but slightly curved along z, as 
displayed at x=4. Due to the strong velocity gra- 
dients in the entrance regions, side layers may 
become unstable exhibiting locally confined 
transient eddy regimes between 0 « x « 1. 


The pressure distribution in the pipe and FCI is 
shown in Figure 4 along the axis (blue) and 
near the sides (red). Pressure gradients in 
MHD flows in bare pipes and in FCI, as pre- 
dicted by an asymptotic theory, are marked as 
straight black lines in the plot. Results from 3D 
numerical simulations perfectly agree with 
these theoretical predictions, the vicinity of the 
FCI entry, center and side pressure separate. 


While the transition of the core pressure ap- 
pears relatively smooth, the pressure near the 
sides experiences a sharp drop in front of the 
FCI. Extremely high current densities entering 
the well conducting wall in front of the FCI re- 
sult in large Lorentz forces, leading to such 
high-gradient peaks near the sides. After en- 
tering the FCI, the side pressure remains al- 
most constant until leveling off with the core 
pressure. For ideal fully developed conditions 
in very long pipes, the relatively low wall con- 
ductance of the FCI in comparison to that of 
the bare pipe wall leads to a significant pres- 
sure gradient reduction factor of about 16. 3D 
MHD phenomena occurring at entrance and 
exit of FCls may reduce the efficiency of the 
insulation due to the additional pressure drop 
Apsp, caused by recirculating 3D currents. 


While experimental results yield data often ex- 
tracted from measurements on the outer sur- 
face of the pipe [3], the presented numerical 
investigations provide valuable data and in- 
sight into 3D MHD phenomena in the fluid. 
Computations confirm the FCI's efficiency and 
show good agreement with asymptotic theory 
and experimental results. 
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Figure 4 Pressure along the pipe axis (blue) and at the 
sides (red). 


MHD velocity distribution and 
pressure drop in manifolds of a WCLL 
TBM 


In the frame of the EUROfusion breeding blan- 
ket research activities, the water-cooled lead 
lithium (WCLL) blanket has been selected as 
"driver" for a DEMO reactor [4]. Test blanket 
modules (TBMs) derived from this concept will 
be tested in the International Thermonuclear 
Experimental Reactor ITER. 


In the present study, liquid metal MHD flows in 
manifolds of a WCLL TBM have been investi- 
gated numerically for different sets of flow pa- 
rameters. 


A view of the WCLL TBM manifold design is 
shown in Figure 5a. It consists of two long po- 
loidal ducts, electrically connected across a 
common wall. The liquid metal flows along a 
series of expansions and contractions due to 
the presence of horizontal stiffening plates that 
separate the breeder chambers arranged in a 
column. This type of MHD flow causes addi- 
tional pressure drop compared to flow in 
straight ducts [5]. Liquid metal flows are inves- 
tigated in the model geometry depicted in Fig- 
ure 5b, whose dimensions are taken from the 
most recent WCLL TBM design. A uniform 
magnetic field is imposed in toroidal direction. 


Numerical simulations are performed for peri- 
odic fractions of the manifolds by assuming a 
stepwise decrease/increase of the flow in the 


feeding/draining manifolds. Velocity and pres- 
sure distributions in the entire manifolds that 
feed/drain 8 breeder zones are then recon- 
structed. For all investigated parameters a uni- 
form flow partitioning among BUs (Breeding 
Unit) has been assumed and flow rates in the 
manifolds have been derived by taking the 
same given mean velocity uo in all breeding 
zones. The pressure drop depends linearly on 
the mean velocity and quadratically on the 
magnetic field strength. Therefore, results can 
be extrapolated and scaled, e.g. to other flow 
rates, depending on design specifications. For 
the discussion of the results, we consider the 
manifold as subdivided into eight parts, each 
corresponding to a specific fed/drained breed- 
ing zone (BZ). For that reason the different por- 
tions of the manifolds along the poloidal direc- 
tion are referred to as BZ1,..., BZ8, as 
indicated in Figure 5a. 


Calculations have been carried out by using a 
finite volume code. Equations describing the 
MHD flow are implemented in the open source 
software OpenFoam [6]. Accurate simulations 
of MHD flows are pretty demanding since they 
require a proper resolution of thin boundary 
layers that form along all walls. Their thickness 
reduces by increasing the intensity of the mag- 
netic field [7]. Moreover, electrically conducting 
walls provide closing paths for electric currents 
and their resolution is crucial to determine the 
total current density in the fluid. About 4-106 
nodes are needed in the two fluid domains and 
6-7-108 in the thick wall to resolve an eighth 
fraction of the manifold. An example of the 
computational mesh used for the simulations is 
depicted in Figure 5c in a toroidal poloidal 
plane. Simulations have been performed on 
the HPC MARCONI! in Italy by using up to 128 
CPUs. 
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Figure 5 (a) Design of WCLL TBM manifold. Feeding and 
collecting manifolds are referred to as duct 2 and duct 1, 
respectively. (b) Periodic fraction of model geometry for 
numerical analysis. (c) Computational mesh for fluid and 
wall. 


Figure 6 (top) shows, as examples for the flow 
at Ha-4000 and uo-0.1 mm/s, 3D views of the 
velocity distribution at three poloidal locations 
in the considered geometry at the bottom of the 
module in BZ1 (a) and in the middle in BZ5 (b). 
In BZ1 the fluid is supplied to duct 2 (feeding 
manifold), while the weak flow in duct 1 (col- 
lecting manifold) is only driven by electromag- 
netic coupling caused by leakage currents. In 
the lower part of Figure 6 the vertical compo- 
nent of the velocity is plotted along the radial 
white dashed lines in the two channels, as 
marked in the 3D views. Due to electromag- 
netic coupling across the common dividing wall 
the flow in duct 2 pulls that in the core of duct 
1 in the same direction. This leads to a buildup 
of pressure along the flow direction that drives 
the backward oriented jets in duct 1 (orange 
curve in Figure 6a at the bottom). In Figure 6b 
results in the middle portion of the manifolds, 
BZ5, are displayed. Here the mean velocities 
in the two coupled ducts are very similar. As a 
result the velocity distribution is analogous in 


both channels. Electric currents induced in the 
cores of the manifolds enter the corresponding 
side walls, they flow towards the external and 
common Hartmann walls through which they 
pass into the opposite parallel wall. Conse- 
quently, no current is exchanged between the 
two ducts and no significant multichannel ef- 
fects are observable. 

(a) B vta © 


0.0435 
0.0365 


B v [m/s] 
0.014 
0.0107 


0.0243 0.0071 
0.0122 0.00356 


-0.0051 0 


Figure 6 Results for Ha=4000 and u=0.1mm/s, at the 
bottom of the module in BZ1 (a), in the middle in BZ5 (b). 
On the upper half of the figure, 3D views of velocity distri- 
bution are displayed at three poloidal locations. On the 
lower half, profiles of vertical velocity are plotted along 
the radial x coordinate, as indicated by the white dashed 
lines. 


In Figure 7a the pressure is plotted as a func- 
tion of the vertical poloidal direction y along the 
central line of the ducts forming the manifolds 
for the flow at Ha=4000 and uo-0.1 mm/s. The 
distance between the two curves, marked by 
red arrows, indicates the pressure difference 
between feeding and draining manifold, i.e. the 
available pressure heads between entrance 
and exit of the corresponding breeder zones. 
This pressure difference is required to over- 
come the pressure drop across the windows in 
the back plate (BP) that connect manifold and 
BUs, and the one along the breeder zone. Ex- 
periments performed to investigate MHD flows 
in a scaled mock-up of a HCLL blanket showed 
that the pressure drop in the BUs, is almost 
negligible compared to the other contributions 
in the blanket module [8]. The additional pres- 
sure drops present at the connections between 
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two portions of the manifolds, i.e. the steps vis- 
ible in the pressure distribution in Figure 7a, re- 
sult from the local 3D MHD effects that occur 
when the liquid metal expands/contracts along 
magnetic field lines. For the flow at Ha-4000 
and uo =0.1 mm/s the total pressure drop in the 
feeding (FM) and collecting (CM) manifolds is 
given by Apru=0.092 bar, Apcw-0.048 bar, re- 
spectively. The pressure drop in the draining 
manifold is smaller than the one in the feeding 
channel because its cross-section is larger. 


pressure head for BU 


T T T T T T T pol tor separating wall 


I 


Figure 7 (a) Pressure plotted along the vertical poloidal 
direction in feeding and draining manifolds for Ha=4000 
and u5-0.1 mm/s. (b) Schematic ideal linear profile of po- 
loidal distribution of pressure in the manifold. (c) Sketch 
of PbLi path and proposed modification of the manifold 
design. 


In an ideal situation, the pressure heads for 
each breeding zone should be equal in order 
to guarantee a uniform flow partitioning among 
all BUs. However, this would be possible only 
if the pressure along both manifolds varied lin- 
early and with same slope, as sketched in Fig- 
ure 7b. In the present manifold design, the 
cross-sections of the two ducts differ from each 
other and remain constant along the poloidal 
length of the module. Consequently, the mean 
velocity in the feeding manifold becomes 
smaller along the poloidal direction, when the 
liquid metal flows from BZ1 to BZ8. Therefore, 
the pressure gradient decreases. Instead, in 
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the collecting manifold the flow rate increases 
along the liquid metal flow path as well as the 
pressure gradient. The larger pressure head at 
the external BUs suggests that the flow rate in 
central BUs is considerably reduced compared 
to the one at the ends of the module. 


Such behavior has been studied in the past by 
experiments performed in a mock-up of a 
HCLL TBM [9] and design modifications had 
been proposed to overcome this undesirable 
effect. A more uniform partitioning of the liquid 
metal among the breeding zones could be 
achieved by adapting the manifold cross-sec- 
tions to the changing flow rates, so that a con- 
stant average velocity is maintained along the 
liquid metal path, and same pressure head 
could be reached for each BU. In order to real- 
ize these conditions, the vertical wall separat- 
ing the two manifolds should be moved step- 
wise in toroidal direction, as proposed in Figure 
7c, where the new wall position is drawn in 
green. 


The present numerical study allows for a first 
estimation of pressure drop in the feeding and 
in the draining manifolds of a WCLL TBM. Nu- 
merical results confirm the typical scaling for 
the pressure drop in inertialess MHD flows un- 
der strong magnetic fields, according to which 
pressure losses scale linearly with the mean 
velocity and quadratically with the magnetic 
field strength, Ap-oLB?uo. However, it is still 
necessary to consider that, even if the pres- 
sure drop in the manifolds represents a large 
part of the total pressure drop in the PbLi sys- 
tem, other contributions should be estimated 
as well. For instance pressure drops in the long 
piping system supplying PbLi to the TBM were 
not part of the present investigation. In these 
pipes PbLi flows under the action of a spatially 
varying magnetic field and the velocity is much 
higher than in the BUs. 


The obtained results show that the pressure 
drop in the manifolds is acceptable for ITER 
applications. However, a more important as- 
pect is the non-uniform partitioning of the liquid 
metal flow among different breeding zones due 
to the geometry of the manifolds in which the 
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cross-sections do not adjust to the changing 
flow rates. The proposed modification of the 
manifold design will be evaluated in future 
studies. 


Design and fabrication of a WCLL TBM 
mock-up for MHD experiments 


For MHD experiments in the MEKKA labora- 
tory a mock-up of the ITER WCLL TBM has 
been designed and fabricated (see Figure 8). 
Objective of these experiments is to achieve 
knowledge about MHD pressure drop, when 
the liquid metal flows around obstacles (water 
pipes) and through contractions and expan- 
sions created by stiffening plates in the mani- 
folds and to determine the distribution of PbLi 
flow among breeder units. 


The WCLL mock-up has been 1:2.5 down- 
scaled to fit into the gap of the magnet at the 
MEKKA laboratory. Figure 9 shows a transpar- 
ent CAD view of the mock-up for visualization 
of the PbLi flow paths. 


PbLi outlet flow 


Figure 9 Transparent CAD view with main dimensions (in 
mm) and PbLi flow paths. 


The WCLL mock-up has been simplified to 
meet manufacturing concerns and to reduce 
fabrication costs, while keeping the geometry 
for the liquid metal flow as foreseen in [10]. An 
exploded CAD view in Figure 10 shows the 
major parts of the mock-up. The main body 
consists of a single piece, which forms the 
walls of all breeder units (BUs), including first 
wall, stiffening plates, back plate and walls 
separating the two inlet and outlet manifolds 
for the PbLi flow. Manifolds are closed by two 
back plates. Details concerning the water 
channels inside the blanket walls have been 
omitted since it should not affect the liquid 
metal flow. Cooling pipes are all present in the 
mock-up and simulated by solid dummy ele- 
ments. All parts are made of 1.4571 austenitic 
steel that has good compatibility with the used 
model fluid NaK. 


Side wall rt 
Baffle plates ~ u 


First wall ~ dummies for 


==+— cooling pipes 
iG 


Main body of 
WCLLTBM ———* 
mockup 


Back plates =— 


PbLi inlet pipe. ———-ga 24 


Side wall —— — — «tX 


Figure 10 Exploded view of the WCLL MHD mock-up with 
major components: main body, side walls, back plates as 
cover for the manifolds, and dummy parts that mimic the 
water pipes. 


The liquid metal is fed into the manifold of the 
main body and removed from it through circu- 
lar pipes, which are connected to the existing 
loop in the laboratory. The PbLi is distributed 
and collected by two manifolds into and out of 
the BUs. Each BU is fed through a small win- 
dow in the back wall of the manifolds. The PbLi 
flow is redirected at the first wall and guided 
back into the outlet manifold (see Figure 9). 


Pressure differences will be measured be- 
tween several points of the mock-up using a 
piping system between pressure taps and the 
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pressure transducers which are located out- 
side of the magnet. Positions of pressure taps 
have been chosen to also satisfy the function 
of feeding, draining and venting (see Figure 
12). 


Flows in different parts of the module are elec- 
trically coupled since currents induced in one 
part may leak across electrically conducting 
walls into neighboring channels. Electromag- 
netic flow coupling and flow partitioning among 
BUs is investigated by the electric potential dis- 
tribution on the walls. Experiments using this 
mock-up in the MEKKA laboratory will start in 
the second half of 2021. 


PbLi inlet PbLi inlet windows PbLi outlet PbLi outlet windows 
manifold into BU 


manifold from BU 


Separation wall Stiffening plates in manifold 


H,O dummies Baffle plates Stiffening plates in BUs 


/ ' | illi W ai 


| "Y 


Figure 11 WCLL TBM mock-up: inlet and outlet manifolds 
(top), breeder zone with dummies representing cooling 
pipes, baffle plates and stiffening plates (bottom). 
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manifold taps pipe 


Figure 12 Finalized mock-up. 


Reconstruction of 3D velocity profiles 
from measured surface potential data 


Liquid metal pipe flow entering a strong mag- 
netic field had been investigated experimen- 
tally and theoretically in a previous campaign 
[11] [12] using the MEKKA facility [13]. Experi- 
mental data for electric potential ® on the ex- 
ternal surface of the pipe had been recorded 
and published in [11]. This data displayed in 
Figure 13 has been analyzed and potential val- 
ues @ on the fluid-wall interface have been 


evaluated by solving a sort of inverse problem 
across the thick-walled pipe. Potential at the 
fluid-wall interface is of particular interest since 
it is a result of the interaction of the flow with 
the magnetic field and from its knowledge it is 
possible to derive an approximate reconstruc- 
tion of the velocity distribution in the entire fluid 
volume. 
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Ha=2000, Re=20026 ‘ 


Figure 13 Distribution of magnetic field B(x) projected to 
the pipe surface and contours of nondimensional electric 


potential D measured on the external surface of the 
pipe. The positions of the potential sensors are marked 
by blue dots. 


A first step in the analysis is a decomposition 
of surface data in circumferential Fourier 
modes 


P(R,a,x) = D(a,x)=> F, (x)sin(ka) , 


where the nondimensional radius R denotes 
the ratio of external to internal radius of the 
pipe and « stands for the angle in a polar coor- 
dinate system, measured with respect to the 
magnetic field direction. The F, are obtained 
as best fit of the approach introduced above 
with the experimental data. Results of F, (x) 
displayed in Figure 14 show that the major 
physical information is already contained in the 


first 3 modes while values of mode 5 are al- 
ready negligibly small. 
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Figure 14 Fourier coefficients F;(x) of surface potential 
(a,x) at the external surface of the pipe for Ha=2000, 
Re-20000. 


Assuming a balance between flow-induced 
electric field v x B and potential gradient V à, 
and the fact that the fluid core potential is con- 
stant along magnetic field lines [14], we can ex- 
press the axial component of velocity as 


u=B 9 _ 106 
O.  Bcosa da ` 


when ¢ is known in terms of its Fourier expan- 


sion. Figure 15 shows velocity contours for the 
two Hartmann numbers Ha=2000 and Ha-5485. 
We observe the formation of increased veloci- 
ties near the sides z=+1 and reduced velocity 
in the center, in the region where the magnetic 
field starts rising. 


Figure 15 Contours of axial velocity u(x,z) estimated 
from measured wall potential for Ha=2000, Re=20026 
and Ha=5485, Re=10043. 


The distributions of axial velocity at the sides 
u(x,z=+1) and in the center u(x,z=0) are plotted 
in Figure 16. For x<-5, both the magnetic field 
and induced potentials are pretty small so that 
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data evaluation for velocity becomes difficult 
and the derived data might become unreliable 
(larger scatter of data for x<-5). However, for 
x» —5 the data shows clear trends. Maximum 
values reach up to uma=1.62 and 2.64 at 
x=-1.68 and x=-3.31 for Ha=2000 and Ha=5485, 
respectively. 


* Ha=2000 
e Ha=5485 


N 


0 Ta T 
-10 5 KO O0 5 10 
eS 


Figure 16 Velocities near the sides (filled symbols) 
u(x,z=+1) and in the center (open symbols) u(x,z=0) for 
two Hartmann numbers. 


lt can be observed that the axial position of 
maximum jet velocity near the sides shifts up- 
stream as the magnetic field Bo increases. This 
is most probably caused by the fact that locally 
the magnetic field increases earlier for higher 
Hartmann numbers. In order to support this as- 
sumption we compare the two experiments in 
terms of their local Hartmann numbers 
Ha(x)=Ha B(x)/Bo as shown in Figure 17, where 
the two positions of maximum jet velocity have 
been added to the figure. It can be seen that in 
both cases the position of maximum jet velocity 
occurs at almost the same local Hartmann 
number near Ha(x)=560. Further downstream, 
the velocities near the sides decrease slowly 
and approach that of a uniform core flow with 
u>0.92. 
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Figure 17 Distribution of local Hartmann number Ha(x) 
for the two cases considered. The marked positions x= 
-3.32 and x- -1.68 denote locations of maximum jet ve- 
locity. 


The transverse distribution of velocity is shown 
in Figure 18 for the high Hartmann number 
Ha=5485 and Re=10043. We observe high ve- 
locity magnitudes near the sides at axial posi- 
tions where the magnetic field starts rising with 
a reduced velocity in the center. Further down- 
stream the flow gradually becomes uniform as 
the magnetic field approached its uniform high 
magnitude. 


The advantage of the presented method for 
evaluation of velocity from measured potential 
data is that it yields a global overview of veloc- 
ity in the entire flow field in contrast to invasive 
measurements where a probe is traversed 
along a single line of a cross-section. More de- 
tails and a discussion about the accuracy of the 
procedure can be found in [15]. 
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Figure 18 Transverse profiles of axial velocity at differ- 
ent axial positions for Ha=5485 and Re=10043. 


Experimental study of heat transfer in 
magneto-convective flows related to 
WCLL blankets 


In WCLL blankets, water-cooled pipes are im- 
mersed in the liquid metal to extract the heat 
from the breeding zone. These pipes generate 
large thermal gradients within the liquid metal 
that may lead to quite complex flows governed 
by the combined effect of buoyancy, electro- 
magnetic force, and pressure gradient. Alt- 
hough some characteristic features associated 
with such flow have already been identified by 
numerical studies [16] [17] experimental inves- 
tigations of magneto-convective flow in proto- 
typical geometries are lacking. Yet, additional 
studies are required to improve the under- 
standing of this type of flow and support the 
complex design of a WCLL blanket. To this 
end, a simplified model geometry was defined 
with the aim of providing a first experimental 
database of MHD buoyant flows with heat 
transfer in WCLL-like geometries. 


Thermocouple probe 


Figure 19 Sketch of the test section. 


A sketch of the test section designed to per- 
form magneto-convection experiments in the 
MEKKA laboratory is shown in Figure 19. It 
consists of the 200x200x100 mm rectangular 
cavity made of PEEK plastic through which two 
parallel horizontal pipes are inserted. Both 
pipes are connected to their own temperature- 
controlled water circuits such that they can be 
kept at constant temperatures 7; and 7» to gen- 
erate the horizontal temperature gradient driv- 
ing the flow. The pipes are made of copper and 
they have an inner core to ensure as good iso- 
thermal conditions as possible. They are 
coated with a thin layer (2 uim) of silicon car- 
bide that provides electrical insulation and pre- 
vents corrosion that could occur at higher tem- 
peratures in contact with eutectic gallium- 
indium-tin, which is used as model fluid. The 
entire test section is instrumented with 71 ther- 
mocouples, 416 electrodes and 3 UDV trans- 
ducers to record simultaneously temperatures, 
electric potentials at the walls, and velocities. 
Additionally, flowrates and temperature 
drop/rise along both pipes are also monitored 
to evaluate the amount of heat exchanged be- 
tween the pipes. More details about the test 
section and the instrumentation can be found 
in [18]. 
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Figure 20 Photo of the instrumented test section par- 
tially thermally insulated before insertion in the magnet 
gap of the MEKKA laboratory. 


Experiments were performed for Hartmann 
numbers up to 3000 and Grashof numbers 
ranging from 2.0-10° to 1.0-10%, corresponding 
to temperature differences AT-T»;- T; im- 
posed between the pipes varying between 
1.4°C and 69.5°C. Whenever possible, the 
temperature of the pipes was chosen such that 
average temperature of the liquid metal 
T = (Tə +T,)/2 was close to the temperature 
inside the magnet gap — set to 30 °C — in an 
effort to limit as much as possible the external 
heat exchanges. 


The temperature distributions were measured 
at the center of the box, on the middle of the 
sidewall (on the heated side, z* = -2.00) and on 
the back plate (x* =-1.98), as well as on the 
middle of the top cover (y* = 1.00) and in both 
pipes. Results obtained are scaled and plot as 
nondimensional temperature 7* defined as fol- 
low: 


, T-T 

T'- Ara 

For hydrodynamic experiments, performed in 
the absence of magnetic field (Ha = 0), a con- 
vection cell forms in the center of the cavity, 
between the two pipes. The buoyant flow re- 
sults in a thermal stratification with the hot fluid 
staying on the top and the cold fluid on the bot- 
tom as seen on the non-dimensional tempera- 
ture profile measured at the center probe 
shown in Figure 21. 
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Figure 21 Nondimensional temperature profiles meas- 
ured at the center of the test section (x* = -0.13 & 
z* = 0) for Ha - 0. 


Magnetohydrodynamic experiments exhibit 
quite different temperature distributions de- 
pending on the strength of the applied mag- 
netic field. At the lowest Ha numbers, the tur- 
bulence suppression by the magnetic field first 
results in larger temperature gradient within 
the flow (Figure 22), i.e. the thermal stratifica- 
tion becomes stronger. As the magnetic field 
further increases, the convection motion is 
damped by the Lorentz forces induced in the 
liquid metal. Consequently, the vertical thermal 
gradient at the center of the cavity progres- 
sively diminishes until the horizontal thermal 
stratification has completely disappeared. 
Eventually, convection motion is so damped 
that thermal conduction emerges as the domi- 
nating heat transfer mechanism and the tem- 
perature distribution becomes symmetric with 
respect to the vertical middle plane of the test 
section at the highest Ha numbers. 


This result is confirmed by the measurement of 
the heat transfer coefficient. The intensity of 
the convective heat transfer has been quanti- 
fied by deriving the amount of heat re- 
leased/absorbed at the hot/cold walls, respec- 
tively, from the energy balance in both water 
pipes. The heat fluxes at the pipes can be then 
evaluated knowing the exchange area and ex- 
pressed in nondimensional form as Nusselt 
number 


rhe, ATL 
u = ——————— 
ndlk(T,; — T) 


Group: Magnetohydrodynamic 


Here, m denotes the mass flow rate of water in 
the pipe and AT, the temperature drop/rise. k 
and c, are the thermal conductivity of galinstan 
and the specific heat of water, respectively, 
and d and / are the diameter and length of the 
pipes (see Figure 23). The characteristic 
length of the problem L considered here is half 
distance between the pipes which is equal to 
the half size of the box measured along mag- 
netic field lines (Hartmann length). 


Figure 23 Sketch for the definition of parameters used 
to derive the Nusselt number. 


The results are presented in Figure 24 where 
Nu is plotted against the quantity Gr/Ha? repre- 
senting the ratio of the buoyancy force over the 
electromagnetic force. The competition be- 
tween buoyancy forces, promoting the convec- 
tion, and the electromagnetic forces, suppress- 
ing the flow, lead to 2 distinct heat transfer 
regimes: 


- For Gr/Ha? < 10, Nu is constant and close 
to unity indicating that heat transfer is 
purely conductive. Almost all convection 
motion has been suppressed by the Lo- 
rentz force. 
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- For 50<Gr/Ha?< 4000, the intensity of 
heat transfer increases with Gr/Ha?. In this 
regime, a scaling law has been established 
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Figure 24 Nusselt numbers measured for several 
Grashof and Hartmann numbers. 


As Gr/Ha? grows further, the pace at which Nu 
increases slows down and eventually the flow 
becomes turbulent as proven by current exper- 
imental results. While this work already pro- 
vides valuable data, extended investigation will 
be carried out focusing on the onset of turbu- 
lence in magneto-convective flow. 


Further work 


In addition to the topics described above, the 
MHD group at ITES KIT contributed in the re- 
porting period 2020 with scientific papers to the 
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Figure 22 Nondimensional temperature profiles measured for Gr- 2-10” and Ha ranging from 0 to 3000 (a) on the 
middle of the top wall (x* = -0.13 & y* = 1) and (b) at the center of the test section (x* = -0.13 & z* = 0). 
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development of fusion technology [19], [20], 
[21], [22], [23], [24], [25]. 
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SAR/ITES research activities in 2020 


Xioayang Gaus-Liu 


In 2020, the research activities of ITES/SAR 
were focused on the analyses and mitigation of 
severe accident of boiling water reactor (BWH) 
and the construction of a new facility ESFR- 
LIVE as a task in the EU Horizon2020 project 
ESFR-SMART. 


In-vessel melt and debris retention is also a de- 
fense-in-depth strategy during the course of a 
severe accident scenario of BWR. The Fuku- 
shima-Accident reminds the urgent necessity 
of a melt retention concept specialized for the 
BWR reactor design. External cooling by flood- 
ing the reactor cavity, the most studied one in 
the past, could be functional for Pressurized 
Water Reactors (PWR), characterized by a 
small reactor pit and low position of the Reac- 
tor Pressure Vessel (RPV) in the cavity. It is not 
suitable for BWR with a deep und large cavity. 
On the other hand, the lower plenum of a BWR 
has a dense array of control rod guide tubes 
(CRGTs), which penetrate the RPV bottom 
vertically to the cavity. CRGTs can contribute 
to the cooling of a relocated debris bed or a 
molten pool in the lower plenum in two ways: 
a) cooling inside the corium by the water- 
cooled CRGTs and b) directing water via 
CRGTs to the top surface of the corium to en- 
able top cooling. The cooling effects of CRGTs 
both in a non-molten debris bed and in a mol- 
ten pool were experimentally simulated in the 
LIVESD test facility [1]. The test configurations 
simulated different scenarios: no cooling at all, 
cooling only by CRGTs and cooling by CRGTs 
together with external wall cooling. 16 CRGTs 
are simulated and were arranged in 4 groups. 
The tubes were mounted on an upper lid with 
water distribution system, and inserted verti- 
cally approaching the bottom of the vessel. 
Figure 1 shows the tube inside.the vessel and 
the water distribution system at the upper lid. 


Figure 1. LIVE-BWR tests with simulants of CRGTs. 


The heating and cooling of debris were per- 
formed in 6 phases, in each of which there was 
a high power period and a low power period, 
as shown by Ptotal in Figure 4, right diagram. 
The external water cooling were only per- 
formed in the phase 5 and phase 6 with the 
combination of CRGT cooling. During the high 
power period, the CRGT cooling could only re- 
move a part of the power, with the conse- 
quence that the debris temperature arise con- 
tinuously. During the low power period, the 
CRGT tubes could effectively reduce the de- 
bris bed temperature and mitigate the thermal 
load on the vessel wall. The cooling of CRGT 
tubes was more effective in a molten pool, as 
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shown in the Figure 2, bottom. The heat re- 
moval rate via CRGT is as effective as the ex- 
ternal cooling. 
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Figure 2. Power input (Ptotal) and heat transfer via 
CRGT tubes (Qrod) and via external cooling (Qwall). 
Top: Debris test, Bottom: molten pool test. 


Another experiment was the investigation in 
melt-down of a debris bed with heterogeneous 
materials. This experiment was assigned by 
the JAEA /CLAD laboratory with the aim to un- 
derstand the accident scenario of Unit 2 Fuku- 
shima accident and thus accomplish safe de- 
commissioning of the Fukushima Daiichi 
Nuclear Power Station [2]. Till present, few ex- 
periments investigate a meltdown process of 
debris bed consisting of materials with different 
melting temperatures. Two experiments were 
performed on the LIVESD facility: LIVE-J1 and 
LIVE-J2 tests. LVE J1 investigates the melt 
down of debris, whereas LIVE-J2 investigates 
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a corium of debris and liquid melt with and with- 
out external cooling. 


The important outcome of LIVE-J1 tests are a) 
the ongoing shrinkage of debris bed volume 
and height during the melt-down of the parti- 
cles with low melting point; b) the shift of 
hotspot in the debris bed and thermal load on 
the vessel wall from bottom toward the upper 
region during the melting. Figure 3 shows the 
debris bed before and after meltdown process. 


Figure 3. Appearance of initial debris bed (top) and the 
debris bed after the test (bottom) in LIVE-J1 test. 


The third accomplishment of SAR group in 
2020 is the design and construction of a new 
large-scale 3D vessel, ESFR-LIVE imitating 
the geometry of the core-catcher of a sodium 
cooled fast reactor (SFR) in a diameter scale 
of 1:6 [3]. The core catcher has the function to 
collect the relocated melt derived from the up- 
per core area and to prevent local thermal at- 
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tack of melt on the reactor vessel. The thermo- 
dynamics of the melt with decay heat and its 
heat transfer at the core catcher boundaries 
will be studied experimentally at KIT and nu- 
merically at CEA. The facility will be in commis- 
sioning in middle of 2021. The ESFR-LIVE test 
vessel simulates the typical tray geometry of 
an in-vessel SFR core catcher, whose cavity 
has a geometry of truncated cone in the lower 
part, and cylinder in the upper part. The inner 
diameter of the upper cylinder is 1 m. 


Fig.4. Top: LIVE experimental platform with ESFR- 
LIVE facility. 1- ESFR test vessel, 2 - heating furnace 
for melt preparation, 3 - pouring spout, 4 - instrumen- 
tation board. Bottom: ESFR-LIVE test vessel with all 
boundary cooling system. 


Four planes of resistance heating element are 
designed to simulate the decay heat of the co- 
rium. Two planes are located in the truncated 
cone (HE1 and HE2) and the other two are lo- 
cated in the cylindrical part (HE3 and HE4). 
The upper two heating planes can be removed 
individually to enable different pool heights at 
the top boundary. The total heating power can 


reach 86.5 kW. The flexibility of the heating 
planes enable a large variety of the pool geom- 
etries, e.g. possible by shutdown the one or 
two of the lower heaters. Correspondingly, bot- 
tom crust grows upwards till the level of the 
lowest heater in operation. The high power 
provided also make it possibly for experiments 
with high decay power, such as debris dryout 
experiments. 


Two experimental series in ESFR-SMART pro- 
ject are foreseen in 2021. The experimental 
data will be used for the validation of the nu- 
merical calculations with TriCAD code in CEA 
and for analytical analysis done by KIT and 
CEA. 


To be published: 
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Detailed investigations on flow boiling of water up to the critical 


heat flux 


Stephan Gabriel, Giancarlo Albrecht, Wilson Heiler, Felix Heineken 


Introduction 


The prediction of two-phase flows with heat 
transfer using system codes and CFD methods 
is essential in various power-engineering appli- 
cations. The current projects of the working 
group Multiphase Flow (MPF) deal with topics 
from the field of In Vessel Retention (IVR) by 
External Reactor Vessel Cooling (ERVC) as an 
emergency procedure as well as safety-related 
evaluation of different designs of small modu- 
lar reactors (SMR) with respect to the boiling 
processes under full-scale conditions. 


The research group concentrates on experi- 
mental investigations, in particular on boiling 
under forced convection up to the critical heat 
flux. The activities in 2020 focused on the 
measurement of the critical heat flux at fluctu- 
ating mass flows, which can occur under 
ERVC conditions at the low-pressure COS- 
MOS-L facility. At COSMOS-H , work was con- 
tinued on completion and commissioning of the 
loops and, in addition, an initial project relating 
to small modular reactors (SMR) was started. 
Within the scope of both projects, the meas- 
urement technology was further developed in 
order to be able to record the relevant meas- 
ured variables in as much as detailed as pos- 
sible. 


COSMOS-L 


Within the scope of the BMWi - project KEK- 
SIMA the critical heat flow under periodically 
fluctuating mass flows is being investigated in 
comparison to steady state conditions. The 
aim is a better understanding of the physical 
processes in IVR/EVRC scenarios as they are 


for example conceivable in Gen Ill+ reactors 
(e.g. AP1000). 


water outlet 


pressure sensor 


windows heated length 


power supply 


pressure sensor 


flow straightener 
water inlet 


Figure 1 CAD-model of the new test section with planar 
heater 


In cooperation with our project partner 
KIT/IATF, a new test section with a flat heater 
(heated length 300 mm, width 20 mm) was de- 
veloped (Figure 1) for the low pressure loop. 
The heater is made of stainless steel and ex- 
tensively instrumented for the measurement of 
temperature distribution, heating power, pres- 
sure gradient and mass flow. In addition, sight 
glasses from three sides allow the use of opti- 
cal measurement technology. This test section 
uses temperature-resistant 3D-printed plastic 
parts for the first time in COSMOS-L. A flow 
straightener and several supports for thermo- 
couples at the heater were manufactured. The 
components can be used permanently at tem- 
peratures significantly above 200°C and have 
already proven themselves in the first experi- 
ments. 
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Figure 2 Simplified P&ID of COSMOS-L with new test section bypass 


In addition, COSMOS-L was extended with a 
bypass parallel to the test section (Figure 2). lt 
is controlled by fast-switching proportional 
valves and has an additional flow meter in or- 
der to measure the mass flow through the test 
section accurately. The mass flow can now be 


Parameter Value 
Pressure 1200 
mbar 
Tinet 60°C 
Subcooling 45K 
Mass flow den- | 185 
sity kg/m?s 
Heating power | 9,3 kW 


varied via the valves without having to change 
the total mass flow of the pump. The profile of 
the variation is thereby programmable. Test 
track and water loop were put into operation 
and a first series of tests was carried out. Fig- 
ure 3 shows two snapshots of the overflowed 


Parameter Value 
Pressure 1200 
mbar 
Tintet 60°C 
Subcooling 45K 
Mass flow den- | 185 
sity kg/m?s 
Heating power | 0 kW 


Figure 3 Instant image of the heated section shortly before the boiling crisis (I) and immediately after the power 


discharge (r). 
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heater from these tests. Both images show the 
side view on the left and the front view on the 
right of the test track respectively. 


First, a number of tests were carried out to de- 
termine the measurement uncertainty of the 
critical heat flux. Subsequently, the first param- 
eter variations were conducted. As the main 
parameter of the investigation, both the period 
and amplitude of the oscillation were varied. 
Experiments were started in which the mass 
flow oscillates periodically (Figure 4). 
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Figure 4 Oscillating mass flow during the first experi- 
ments 


During the tests, the temperature distribution 
was observed through thermocouples at- 
tached to the back of the heater. The boiling 
crisis was detected by the associated temper- 
ature excursion at the heater and the heating 
power was dropped to protect the test section 
(Figure 5). 
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Figure 5 Temperature distribution at the heater 


The values obtained this way already show a 
measurable influence of the fluctuations in the 
first tests. In Table 1 the example of two se- 
lected measuring points shows a reduction of 
the critical power with fluctuating mass flows of 
approx. 2.7% with a measurement uncertainty 
of approx. 1%. This already occurs with a mass 
flow amplitude of 30% of the nominal value; in 
the postulated scenario of an IVR, much larger 
amplitude values can also occur. 


Several test series with additional optical in- 
strumentation and parameter variations will be 
conducted with this test setup over the next few 
months. In order to investigate larger mass 
flows with constant inlet subcooling, the ther- 
mal power of COSMOS-L facility will be in- 
creased in the near future by an additional heat 
exchanger in order to be able to investigate 
larger mass flows with constant inlet subcool- 
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Figure 6 Extract from the test evaluation of the first test series. Comparison of a measurement with stationary mass 


flow and a measurement with fluctuating mass flow 
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COSMOS-H 


Complementary to the activities on the low- 
pressure facility COSMOS-L, the completion 
and commissioning of the high-pressure facility 
COSMOS-H is currently in progress. In 2020, 
further subsystems of the plant were finalized. 
Figure 7 shows the simplified flow diagram of 
the plant, which has two cooling circuits in ad- 
dition to the high-pressure circuit. With an in- 
stalled thermal power of approx. 2 MW and a 
working pressure of up to 17 MPa at approx. 
360?C, the water/steam circuit will reach power 
plant conditions in test operation. 


In the EU project McSafer, which started in 
2020, investigations into the safety of various 
SMRs [1] are planned together with the project 
partner KIT/INR and other partners from re- 
search and industry. In COSMOS-H, experi- 
ments on heat transfer from nucleate boiling up 
to the critical heat flow will be carried out on 
two different test setups, a single tube in the 
annular gap and a rod bundle consisting of five 
zircalloy-4 tubes. Comparable experiments 
have already been performed with a smaller 
heated length under low pressure conditions at 
COSMOS. [2, 3]. 


Figure 8 shows the first test arrangement. The 
heated section was constructed and the first 
components have already been manufactured. 
The assembly will be fully built up in the follow- 
ing months. In addition to thermocouples and 
pressure sensors, laser-based glass fibre sen- 
sors [4] will also be used in these experiments 
to measure the phase distribution. 


To create a well-defined inflow condition, a 
flow straightener will be positioned at the inlet 
of the annular gap. Caused by the electrical, 
thermal and hydraulic requirements, this part 
was realized as a ceramic 3D print [Figure 9]. 


Furthermore, the test track has several sight 
glasses so that optical methods such as high- 
speed videometry can also be used in future 
experiments. As the test track has an in- 
creased risk of leaks due to the glass windows, 
a multi-stage safety concept was created for 
the plant. For this reason, a fragment shield 
was constructed around the test section to 
catch glass splinters and the steam jet in the 
case of a window failure [Figure 10]. 


Figure 7 Simplified P&ID of COSMOS-H 
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Figure 8 Test arrangement with a single tube in an annular gap 
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Figure 9 Flow straightener in CAD model and as 3D printed part made of AL2O3 ceramics 


Figure 10 Pressure hull of the test section including high-pressure windows (I) and fragment shields (r) 
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Introduction 


Research and development activities at the 
Karlsruhe Liquid Metal Laboratory (KALLA) in 
2020 can be categorized in three pillars: 


- Thermal fluid dynamics of liquid metals 
- Liquid metal based heat storage 
- Liquid metal based process technology 


In all three areas, experimental and modelling 
approaches are combined, mutually benefitting 
from each other. All three areas do also benefit 
from the long lasting liquid metal engineering 
experience at KALLA and contribute to its pro- 
gress at the same time. This holds for the de- 
sign of the loop facilities, reactor setups, instru- 
mentation, components as well as for process 
control and data acquisition systems. They all 
define the state of the art in liquid metal tech- 
nology and enable innovative applications. The 
latter range from liquid metal based high tem- 
perature thermal storage up to highly innova- 
tive process routes for e.g. the production of 
energy carriers. 


High temperature heat storage is expected to 
be an enabling technology for post-fossil heat 
supply to high temperature thermal process 
technologies, based on fluctuating sources or 
waste heat. In terms of high temperature pro- 
cesses itself, liquid metals become more and 
more visible as heat transfer or even reaction 
media in highly innovative process routes like 
direct dehydrogenation of methanol and direct 
thermal decomposition of methane. While the 


first can be part of a process chain to produce 
sustainable liquid energy carriers, the second 
can be used as hydrogen production technol- 
ogy or as part of innovative pathways towards 
utilizing carbon dioxide as a raw material for 
the production of e.g. carbon black. All of the 
mentioned topics fit together in a picture of en- 
abling technologies and new ideas for the en- 
ergy transition with its challenges not only for 
the energy sector, but for chemistry, transport 
and many others. 


In terms of serving fundamental research, a re- 
furbishment of the GALINKA-loop has been 
accomplished, that will facilitate precise meas- 
urements of the heat transfer behaviour in non- 
homogeneously heated tubes, delivering both 
— correlations as well as validation data for so- 
phisticated simulation studies and turbulence 
model development. A second important ex- 
perimental topic with direct relevance for simu- 
lation tool development, validation and safety 
assessment of innovative nuclear research in- 
stallations is the investigation of local block- 
ages and deformations in liquid metal cooled 
rod bundles. Here we continue the traditional 
combination of generic investigations with ap- 
plication relevant parameter ranges, resulting 
in unique data sets that are acknowledged 
worldwide in nuclear safety research. 


We are grateful that - despite the difficulties of 
the pandemic situation - we could continue the 
strong partnerships with researchers across 
Europe and worldwide, and even establish 
new partnerships with internationally leading 
companies in the fields of liquid metal and pro- 
cess technology. Please see brief descriptions 
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of all mentioned activities on the following 
pages and do not hesitate to contact us, if you 
see an opportunity to develop creative ideas 
with the KALLA team in future. 


Thermal Fluid Dynamics of Liquid 
Metals 


GALINKA — Gallium-Indium-Tin loop 
Karlsruhe 


One of the most traditional installations in 
KALLA, the liquid metal test loop GALINKA, 
has been completely renewed in order to real- 
ize high-precision liquid metal thermohydraulic 
experiments at room temperature. The used 
eutectic composition of gallium, indium and tin 
has a melting point of 10.5?C. The loop is 
shown on the left-hand side of figure 1. It hosts 
a variety of instrumentation to determine mass 
and volume flow rate, temperatures and pres- 
sures at several places in the loop. Experi- 
ments with heat loads up to 20 kW, which can 
be cooled back with a cooling water circuit and 
a liquid metal volume flow rate up to 1.5 m?/h 
can be realized. Due to the low operating tem- 
perature and the low installation space of the 
loop compared to other liquid metal loops, new 
experiments can be realized easily and cost ef- 
fective. 


In a current project supported by the German 
Research Foundation (DFG) (WE 4672/4-1), 
the influence of an asymmetric thermal bound- 
ary condition on the heat transfer in a turbulent 
liquid metal pipe flow is under investigation. 
This type of boundary condition can be found 
in concentrating solar power plants [1, 2]. For 
this purpose, a test section has been devel- 
oped to realize either symmetric or asymmetric 
thermal boundary conditions. Pacio et al. [3] al- 
ready showed, that even for a symmetric ther- 
mal boundary condition high-precision thermo- 
hydraulic experiments would be beneficial. 
Prior to the liquid metal experiments, the new 
test section has been validated using water at 
different Prandtl-Numbers and the results were 
compared to Gnielinski's correlation for heat 
transfer in turbulent pipe flows [4] showing a 
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good conformity. An according parity plot can 
be seen figure 1. 
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Figure 1: GALINKA loop with all essential components 
(top), validation of testsection using water experiments 
and Gnielinski's correlation for heat transfer in turbu- 
lent pipe flows (bottom). 
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Fuel assembly experiments in Heavy liquid 
metal flows 


For the safe and reliable operation of a fuel as- 
sembly in the reactor core, the knowledge of 
the heat transfer to the coolant is essential. 
Moreover, during the lifecycle of the assembly 
its geometry can be deformed by swelling, 
creeping and mechanical defects or blocked by 
debris. As a consequence, locally reduced 
cooling and hot spots are expected. 


A set of two different experiments and accom- 
panying CFD simulation is planned in support 
of safety studies for heavy liquid metal cooled 
fast reactor systems in the EU-Projects PATRI- 
CIA and PASCAL. In a first study, the effect of 
a well-defined porous blockage in a wire 
spaced 19-pin rod bundle will be investigated. 
Following the EU-Project MAXSIMA, where a 
total blockage has been investigated, the effect 
of a more realistic sintered blockage with well 
defined porosity will be investigated in PATRI- 
CIA. Detailed instrumentation will give insights 
about local hot spots as well as recirculation 
patterns which will be used for the validation of 
numerical models in house by the Framatome 
Professional School (FPS) Group of the ITES 
and external partners in the framework of the 
EU-project. In a second study, the effect of de- 
formations in a rod bundle on the local flow and 
temperature field will be studied. In the frame- 
work of EU-Project PASCAL, a water rod bun- 
dle experiment will be set up at the Karman In- 
stitute for Fluid Dynamics (VKI), Brussels with 
design and numerical support by KALLA. 


Figure 2: a) Design detail for the planned rod bundle 
deformation experiment in EU-Project PATRICIA b) 
Design detail for the planned rod bundle blockage with 
well-defined porosity in EU-Project PATRICIA 


Figure 3: Complete blockage of two sub channels for the 
experiments in EU-Project MAXSIMA, wires are part of 
the temperature sensor instrumentation within the block- 
ages and the surrounding rod walls. 
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Velocity: Magnitude (m/s) 
1.0072 1.5109 


0.00000 0.50362 


Figure 4: Velocity magnitude plotted on selected planes 
near the blockage and the streamlines colored with veloc- 
ity magnitude for the blockage experiments in 
MAXSIMA!. 
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Liquid Metal Based Heat Storage 


Liquid metals are tested for the use in thermal 
energy storage systems at Karlsruhe Liquid 
Metal Laboratory (KALLA), as they facilitate 
highly efficient heat transfer in a wide temper- 
ature range. A packed-bed storage system 
with filler material was proposed as a promis- 
ing storage system in previous theoretical 
works at KALLA [1,2]. Starting in 2020, experi- 
mental tests are carried out, following a series 
of subsequent experiments with increasing 
complexity, size and technology readiness 
level during the upcoming years. First, the op- 
erability of a liquid metal thermal energy stor- 
age with a packed-bed has been successfully 
shown in the lab-scale experiment VESPA [5] 
with spherical ceramic particles as filler mate- 
rial and lead-bismuth as the heat transfer fluid 
[3,4]. Figure 5 shows a picture of the open stor- 
age container filled with zirconium silicate par- 
ticles. 


The material of the filler material has been se- 
lected in cooperation with the group of Prof. 
Müller and Dr. Weisenburger at the Institute for 
Pulsed Power and Microwave Technology 
(IHM) at KIT by storing various candidates (ce- 
ramics, glasses and natural stones) in lead- 
bismuth at very high temperatures in order to 
examine their corrosion behaviour via optical 
and electron microscopy and other advanced 
characterization methods. Important lessons 
regarding filling, draining, cyclic operation, 
measuring of temperature and flow have been 
learned for the planned next experiment, fea- 
turing a demonstrator with 100 kWh thermal 
energy storage capacity that is to be built and 
tested at KALLA. 


Figure 5: Packed-bed thermocline storage in the pre-ex- 
periment VESPA at KALLA (Photo: F. Müller-Trefzer) 
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Liquid metals are especially suitable for high 
temperature applications, as they are liquid up 
to temperatures of more than 1000°C. How- 
ever, their corrosive nature requires technical 
solutions regarding the structural materials that 
can be used, especially at temperatures be- 
yond 600*C. In this regard, a BMWi-project to- 
gether with the German Aerospace Center 
(DRL) and the industry partner KSB, an inter- 
national producer of pumps and valves, is 
starting in 2021 regarding testing of pumps and 
valves at temperatures up to 700?C [6], again 
also including the colleagues from the group of 
Prof. Müller and Dr. Weisenburger at IHM. 


References: 


[1]  Niedermeier, K., Flesch, J., Marocco, L., 
Wetzel, Th. (2016) Assessment of thermal en- 
ergy storage options in a sodium-based CSP 
plant. Applied Thermal Engineering, 107, 386- 
397. 


[2] Laube, T., Marocco, L., Niedermeier, K., 
Pacio, J., Wetzel, T. (2020) Thermodynamic 
Analysis of High-Temperature Energy Storage 
Concepts Based on Liquid Metal Technology. 
Energy technology, 8 (3) 


[3] Niedermeier, K.; Müller-Trefzer, F.; Wei- 
senburger, A.; Wetzel, T. (2020) High-temper- 
ature storage with liquid metals - Design of a 
prototype storage system and material testing. 
2nd International Workshop on Carnot Batte- 
ries (2020), Online, 15.-16. September 2020 


[4] Niedermeier, K.; Müller-Trefzer, F.; 
Daubner, M.; Marocco, L.; Weisenburger, A.; 
Wetzel, T. (2020). Theoretical and Experi- 
mental Studies of Dual-Media Thermal Energy 
Storage with Liquid Metal. 26th SolarPACES 
(2020), Online, 28. September-2. Oktober 
2020 


[5]  Dieterle, K.; Inbetriebnahme und experi- 
mentelle Charakterisierung eines Wärmespei- 
chers mit Flüssigmetall als Wärmeträgerme- 
dium, M. Sc. Thesis. Karlsruhe Institute of 
Technology (KIT), 2020 


[6] Heidelberger, M. (2021) Erneuerbare 
Energien: Auf dem Weg zu thermischen Groß- 
speichern. KIT press information 034/2021. 


Liquid Metal Based Process 
Technology 


DECAGAS™ - Decarbonisation of nature gas 
in liquid metal 


Hydrogen technology is seen as a key to the 
success of the energy transition. In June 2020, 
the German government adopted the national 
hydrogen strategy [1]. These extensive 
measures and investment packages are in- 
tended to establish hydrogen as a predominant 
energy carrier in the foreseeable future. 


In recent years, intensive work has been car- 
ried out at KALLA on a continuously operable 
process for direct thermal methane pyrolysis 
as a viable process to produce hydrogen and 
solid carbon from methane and - perspectively 
— biogas and natural gas. In the process, the 
extensive experience with liquid metals, estab- 
lished at KALLA, could be transferred into suc- 
cessful lab-scale experiments, which to date 
have repeatedly proven the possibility to reach 
high methane conversion rates and their de- 
pendency on operation parameters like tem- 
perature [2, 3]. 


Together with the industrial partner Wintershall 
Dea [4], the technique is now being further de- 
veloped for the use on an industrial scale. The 
basis is the pyrolysis of methane in a liquid 
metal bubble column reactor (see Fig. 6). In 
this process, liquid tin is flown through by gas- 
eous methane at 800 *C - 1200 *C. Due to the 
temperature conditions, the methane is split 
into two components, solid powdered carbon 
and gaseous hydrogen. Since the solid carbon 
has a much lower density than the liquid tin, it 
accumulates on the surface of the liquid metal 
and can be removed from the reactor. The gas 
phase, containing a mixture of hydrogen and 
unreacted methane, will be further processed 
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according to the requirements of the final ap- 
plication. Unreacted methane and other hydro- 
carbons can be fed back into the reactor. 


© 


Erdgas 
Biogas 


Figure 6: Graphical illustration of methane pyrolysis using 
bubble column reactor. (Infographic: Leon Kühner, KIT) 


[4] 
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NAMOSYN - Sustainable mobility with 
synthetic fuels 


The BMWi project NAMOSYN (Nachhaltige 
Mobilitát mit synthetischen Kraftstoffen) was 
created and started in 2019 [1]. Its main aim is 
to develop and test synthetic fuels that can be 
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produced and used sustainably. One promis- 
ing possibility is to use formaldehyde as an in- 
termediate to synthesise, e.g. OME, an alter- 
native to diesel. Moreover, it is also by far one 
of the most commonly produced chemicals in 
the industry. Nevertheless, a stoichiometric 
amount of water is formed in the current meth- 
ods of CH20 production and needs to be sep- 
arated, which makes these processes costly. 


In contrast, the direct dehydrogenation of 
methanol with sodium vapour as catalyst may 
allow to avoid the energy-consuming water 
separation process and to obtain highly con- 
centrated formaldehyde [2]. Despite these po- 
tential advantages, there is no industrially ma- 
ture method for the production of anhydrous 
formaldehyde yet. The MEDENA set-up is be- 
ing constructed at KALLA in order to prove or 
develop according process concepts [3]. A 
number of parameters, such as the amount of 
sodium, the mole fraction of methanol and tem- 
perature, play a key role in this type of reaction 
and will be investigated in detail. The MEDENA 
set-up aims to demonstrate that vaporised ele- 
mental sodium homogenously catalyses the 
dehydrogenation of methanol to anhydrous for- 
maldehyde with high activity and selectivity. 


Figure 7: MEDENA set-up for the dehydrogenation of 
methanol to anhydrous formaldehyde at KALLA (Photo: 
M. Kamienowska) 


Group: Karlsruhe Liquid metal LAboratory (KALLA) 


References 


[1] X Project website: http://namosyn.de/ 


[2] Ruf, S.; May, A.; and Emig, G.; "Anhy- 
drousformaldehyde by sodium catalysis", Appl. 
Catal. A. Gen., vol. 213, no. 2, pp. 203-215, 
May 2001. 


[3] Cortelezzi, D.; Auslegung einer Ver- 
suchsanlage zur trockenen Dehydrierung von 
Methanol zu Formaldehyd, M. Sc. Thesis. 
Karlsruhe Institute of Technology (KIT), 2020 


NECOC - Negative Emissions by Separation 


of atmospheric CO» into technically usable 
Carbon Black 


The BMWi funded project NECOC (Negative 
Emissions by Separation of atmospheric CO» 
into economically usable Carbon Black) fo- 
cuses on greenhouse gas reduction and circu- 
lar economy based on CO». In collaboration 
with the two start-ups Climeworks Germany 
GmbH and INERATEC GmbH a group of re- 
searchers from the Karlsruhe Liquid Metal La- 
boratory (KALLA) and the Institute of Thermal 
Process Engineering (TVT) is building an inno- 
vative demonstration plant [1]. The objective of 
the project is to combine three process steps 
to transform atmospheric CO» into economi- 
cally usable, solid carbon. 


Atmospheric CO: is seperated from the air via 
a direct air capture (DAC) facility developed by 
Climeworks. The second process step, oper- 
ated by INERATEC, is the catalytic methana- 
tion of CO2 using microstructured reactor con- 
cepts in order to overcome the challenge of the 
high thermal stability of the CO2 molecule. In a 
third step the thermal dissociation of methane, 
also known as methane pyrolysis, takes place 
in a liquid metal bubble column reactor devel- 
opped at KALLA. The innovative technology 
using liquid tin as heat transport medium pre- 
vents the pyrolysis reactor from clogging as the 
solid carbon produced by methane pyrolysis 
rises to the liquid metal surface as a powder 
due to its lower density compared to tin. 


Besides demonstrating the feasability of the 
combined process another objective of the pro- 
ject is a closer analysis of the pyrolytic carbon. 
Therefore, and for the continuous operation of 
the process, a modification of the pyrolysis re- 
actor allowing continuous carbon recovery on 
lab scale has been developed and tested dur- 
ing the last months. At the same time construc- 
tion works for the demonstration plant started 
several weeks ago after detailed planning. In 
addition to the BMWi funding, an equipment 
application was successful at KIT within the 
framework of the Excellence Initiative (KIT fu- 
ture fields) and as a result a very well equipped 
gas chromatograph could be procured. This 
will make it possible to identify by-products in 
the product gas stream across a broad spec- 
trum. 


CO--haltige Luft 


Methanspaltung 


u 


Methanisierung 


CO:-Abscheider 


CO--freie Luft 


Elektrolyse 


vom klimaschädlichen Treibhausgas CO: 
^. zum Hightech-Rohstoff Carbon Black 


=C0.C 


Figure 8: Schematic of the process for converting the 
greenhouse gas CO2 into technically usable carbon 
black, which is currently being set up at KIT in a plant 
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Smart Resilience Engineering: Mitigative Planning of 
Networked Infrastructures, Uncertainties and Robust 


Management 


S.S. Ottenburger, W. Raskob, T. Müller, S. Móhrle 


Introduction 


At present, systemic supply risks in a future 
smart urban world, where networking, automa- 
tion and complexity of infrastructure systems 
have increased, can hardly be quantified relia- 
bly [1]. Part of this uncertainty is that the long- 
term steady planning and reliable operation of 
sustainable energy systems, which depend on 
volatile renewable decentralized feed-in, are 
based, among other things, on regional climate 
or weather forecast models. However, the fur- 
ther we look into the future, the more uncertain 
these are. 


For instance, an additional bulk demand for 
electricity for the use of cooling systems asso- 
ciated with a heatwave [2] can lead to supply 
bottlenecks and overloads of renewable distri- 
bution and transmission grids, ultimately caus- 
ing large-scale blackouts. The resulting failure 
cascades in highly networked and automated 
or electrified systems [3] may lead to consider- 
able supply failures of systemic proportions 
and significantly reduce the performance of 
complex critical infrastructure systems as typi- 
cally existing in urban environments. If resili- 
ence is not integrated into a system design, re- 
newable energy systems run risk to fail in the 
long-term. Thus, sustainability and resilience 
are thought together; this working group devel- 
ops new methodologies and concepts for iden- 
tifying resilient planning patterns and smart 
management systems for adaptive and auton- 
omous energy grids resp., establishes sys- 
temic risk assessment frameworks and early- 


! https://portal.iket.kit.edu/projects/MCDA/ 


warning concepts, and furthermore designs 
decision support systems for smart crisis man- 
agement in complex multi-stakeholder environ- 
ments. 


These research topics belong to the emerging 
field of ‘Smart Resilience Engineering’, which 
leverages lived principles and long-standing 
expertise of this working group: JRODOS is a 
decision support system for nuclear emer- 
gency management, which has been devel- 
oped and hosted by this group for many years 
and is operationalized in more than 40 coun- 
tries worldwide. However, the gain in 
knowledge in nuclear emergency management 
over the last decades is not yet at an end, as 
on a global level, it can be assumed that nu- 
clear power generation is more likely to be ex- 
panded in order to achieve climate targets, and 
therefore, for example, in the case of Small 
Modular Reactors (SMRs), which would also 
be installed close to cities and dense critical in- 
frastructure systems, research on the safety of 
SMRs and emergency protection must be 
adapted to new siting concepts. 


An example of the synergetic effects is the 
MCDA-KIT! [4-6], a general applicable multi 
criteria decision analysis (MCDA) framework, 
has been developed by this group in the con- 
text of framework, that has been nuclear emer- 
gency management (Fig. 1). Originally created 
for JRODOS, the MCDA-KIT applies MCDA to 
support decision making in emergency and 
preparedness, especially if the process is not 
distinct as measurements and simulation are 
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always affected by uncertainty. On top of that, 
as in such situations decision-making is in gen- 
eral also supervised by an advisory body, peo- 
ple with different backgrounds and expertise 
are involved in the process, causing different 
opinions to collide. Finding a robust and ac- 
ceptable consensus for everyone involved is a 
difficult task. In this context, MCDA as a sup- 
porting method has become increasingly pop- 
ular. The MCDA-KIT has successfully been ap- 
plied and enhanced in many European 
projects [4—6], qualifying the group as advisory 
for IAEA. Moreover, further collaborations 
were established, e.g. with the group "Sustain- 
able Biveconomy” at the Institute for Technol- 
ogy Assessment and Systems Analysis at KIT 
developing an over-institutional tool, where the 
main objective is to create a framework for in- 
tegrated sustainability and resilience assess- 
ment. 


JRODOS: Smart Resilience of 
Nuclear 


en Energy Systems, 
Emergency Resilience 


2 Critical 
Management Engineering infrastructures 


Figure 1 Research in Smart Resilience Engineering de- 
velops new methodologies benefiting from long-standing 
expertise in nuclear emergency management, while in 
turn adapting modern concepts of resilience engineering 
leads to innovations in JRODOS; some key competen- 
cies of these domains are highlighted in the surrounding 
boxes. 


In conclusion, by combining design questions 
and management options, this group repre- 
sents challenging research topics in resilience 
engineering that deal with different types of un- 
certainty — thereby, smart methodologies, 
models, and algorithms are developed, also 
applying artificial intelligence, in order to iden- 
tify robust solutions be prepared for unex- 
pected events, to mitigate adverse impacts 
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through resilient designs of networked infra- 
structures, and to support decision making in 
effective crisis management (Fig.1). 


Resilient Smart Grid Design and 
Microgrids 


Basically, we believe that microgrids (MGs) 
contribute to a more resilient energy system in 
view of the aforementioned risks, see e.g. [7]. 
Applying criticality-based decision criteria for 
clustering MGs on the distribution level is a 
promising approach to link critical infrastruc- 
ture protection to the energy system. Integrat- 
ing this principle into network planning appears 
to be a feasible way forward in the context of 
energy system transformation, in which large 
investments are expected anyway [8]: this is 
not about extra measures to strengthen resili- 
ence, but about understanding certain topolog- 
ical degrees of freedom in design to increase 
urban resilience. 


In this context, MGs can be seen as built-in re- 
dundancies giving rise to a cellular clustering 
of a distribution area: local resources are able 
to be self-sufficient in a small area for a certain 
period of time and furthermore MGs can switch 
to island mode if the overall grid fails to supply 
power [9]. How to design MG structures is not 
cast in stone as well [10—12]. 


Roughly speaking, the size of a MG depends 
on the power sources, storages that are in- 
stalled, the power consumption profiles of the 
infrastructures/consumers in there, and the 
maximum expected duration of self-sufficient 
power provision. Thus, topological degrees of 
freedom appear in the design of power and 
smart metering infrastructures but as well in 
terms of cellular clustering of an energy sys- 
tem. 


There are known graph-based network metrics 
that can be applied as weighted indicators in 
composite resilience metrics by making use of 
multi-criteria decision analysis (MCDA), where 
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Table I: A list of known graph-based network metrics see e.g. [13] 


EY characteristic path-length Efficiency- average of the shortest path lengths 
EN central point dominance Dominance of particular nodes 


f critical ratio of defragmenta- 
tion 


Ch 
c 


algebraic connectivity 


Removal 
defragmentation of the network in different 
clusters — robustness of a network against cat- 
astrophic natural events (e.g. earthquake) 


of a fraction 2f. leads to a 


Expresses the number of disjoint paths, i.e. the 


network remains fully connected despite the 
removal of nodes 


meshdness coefficient 


Robustness and path redundancy - failure of 


nodes 


the weights can be specified through an ana- 
lytical hierarchical process (AHP) with the help 
of stakeholder engagement [13]. 


If we consider MG clustering, statistical peak 
load values and initial criticality of network 
nodes as further attributes of a distribution net- 
work, there is the following metric analogously 
defined as in [14]: 


Let M be the set of all MGs in a smart grid (SG) 
and J, be the index set of all infrastructures 
that lie in an MG A € M. Furthermore, we de- 
note the statistical power peak demand of a 
consumer i € J with d;. A measure for critical- 
ity-demand concentration is given by 


ED” a NUN DE 
A = Xjes (5 =) "(zz x) e 


where 0 € x,y € 1 and x + y = 1. Depending 
on the values of x and y, the greater this indi- 
cator is the more critical and tense the supply 
of infrastructures in A and the provision with 
critical services might get in times of disturbed 
power supply. The new metric associated to a 
SG is given as 


CD :— max CD, , (2) 


where x and y where omitted for better reada- 
bility. 


Hence, CD can be integrated into an MCDA- 
based resilience metric as described above for 
urban resilience assessments of smart distri- 
bution grids. It is important to point out, that the 
list of metrics given here does not refer to 
power physical constraints, however they play 
a crucial role from a network perspective. Nev- 
ertheless, power physical constraints or met- 
rics should be considered to better assess the 
feasibility of power system designs - this is 
subject to current research. Furthermore, this 
approach can be adapted to any other net- 
worked infrastructure including hydrogen sup- 
ply systems. 


Resilient Energy Management 


The IEEE standard 1366 [15] specifies a series 
of twelve reliability indices for power distribu- 
tion systems, the most familiar of which are 
SAIDI (System Average Interruption Duration 
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Index) and SAIFI (System Average Interruption 
Frequency Index). Roughly spoken, SAIDI and 
SAIFI refer to periods without power supply 
and are useful in a posteriori assessments of 
the reliability of power distribution systems. 
LOLE (Loss of Load Expectation) refers to the 
number of hours per year during which supply 
is statistically expected not to meet demand in 
the long run. Since these indices are deter- 
mined statistically, they are not subtle enough 
for the objective of having fine resolved meas- 
urements of security of supply as a basis for 
short-term or real-time decision-making. 


The purpose of this subsection is to propose a 
(global) criticality-based metric that might be 
used in some sort of composite resilience met- 
rics also taking efficiency and fairness into ac- 
count [16]. 


We assume that infrastructures are equipped 
with smart meter, which are able to prevent 
power consumption above a specific threshold 
and furthermore have a certain power demand 
flexibility that depends on their sub-processes 
and the current demand of the infrastructure's 
functions. If J := {1,...,N} is the set of all infra- 
structures from an urban distribution network, 
then we express the demand flexibility of an in- 
frastructure lel with the interval 
[P$ min Ph max] and its global criticality with cı. 
Of course, there are infrastructures that need 
to maintain all sub-processes in order to fulfill 
their main tasks sufficiently, e.g. in dialysis clin- 
ics, dialysis machines, pumps, hot water disin- 
fection systems etc. have to be in operation 
mode concurrently. 


Let us consider k intervals of fixed length, e.g. 
30 minutes, and let sp, , be the assigned power 
to infrastructure l in time interval t, which is 
considered as the power threshold l cannot ex- 
ceed. We define the following urban resilience 
indicator for power distribution: 


si := Xin Xa & * Qi, (3) 
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where — assuming there is at least one i € 


(1, ..., N} with c; #0 such that &:= z4"— is 
XMj-i6j 


well-defined — c; is called the normalised vari- 
ant of c; and 


l 
0 , Spit < Pb min 
SPLt l l 
pl , Po min < Spit < Pp max (4) 
D,max 


Que: = 
l 
1, Pp max < Spit 


characterizes the quality of power supply from 
the customer's perspective. In [16] first simula- 
tion results were presented, which have shown 
that instead of having controlled load shedding 
as described in Fig. 3, operationalizing this 
type of indicator in SGs can result in urban re- 
silient — fair and efficient — power flows avoid- 
ing large-scale blackouts and in a significantly 
improved si, see Fig. 2. 


Instead of having a controlled load shedding as 
described in Fig. 3, where large areas were 
switched off in a round-robin manner, smarter 
distribution could avoid such blackouts and as- 
sign power thresholds in an urban resilient and 
physically feasible way to the customers - i.e. 
no large-scale blackouts would occur, see Fig. 
2. 
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Figure 2 Heat map of the relative power supplies: in contrast to the rolling blackout of Fig. 3 continuous power flows 
are realised in an urban resilient, fair and efficient way over four time steps (red colour indicates induced blackout) 


[16] 


"n 
Hiii 


in HHH 
— HW 


m 
I HHHH 
"rro 

— HH 
r1 
iiia didi 
Mrrnratraaiaainna 


— HH 


"n 
$ 
= 

E 
o 
a 
d 
m 
<= 
to 

c 
2 
z 

© 

n 
a2 

o 
E 

v 
E 
5 

o 
8 

ap 
A 

= 
2 

D 

e 


a 
Hii 


IHH HHH 
Ln 
m 
IHH 
rennen HHH 
Hii 
| rtt 
0 -HHHH 
GrrrirÜ ir it iil 


His 


pertum 
| iH 


AH 
HHH 
¡EPT ag 
iii 
HHH 
egy naa 
iii 


1 z 3 1 3 i 1 1 


Figure 3 Three supply areas alternately taken off the grid in a round-robin manner — red colour indicates induced 


blackout [16] 
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Dynamic Measure Development — 
Agent-based optimization 


Well-known features of agent-based modelling 
include scalability and simple adaptability of 
specific agent models and parameters. Be- 
sides a flexible integration of new agents, also 
further preferences can be included on an 
agent-specific level and calibrated or dynami- 
cally weighted depending on the agent's envi- 
ronment and forecasts w.r.t. various aspects 
related to the demands and resource availabil- 
ity e.g. weather forecasts in conjunction with 
renewable energy generation. Furthermore, 
based on the various needs of agents, repre- 
senting different stakeholders or entities in a 
system, their preferences, and global strategic 
goals, solutions on how to distribute resources 
in a resilient way can be decided by negotia- 
tion- or optimization procedures targeting at 
identifying feasible and robust measures, and 
simulated. Agents' preferences and global 
strategic goals are considered as a result of a 
two-stage MCDA-process, which might be in- 
fluenced by boundary conditions or sudden 
state changes. 


Therefore, a resource independent agent- 
based optimization framework (ABO-F) has 
been set up that distinguishes between three 
different agent types: networked-, non-net- 
worked-, and purely demanding agents. Once 
this ABO-F has been instantiated according to 
a specific context, we speak of a multi-agent 
system (MAS). This MAS is equipped with a 
certain resource taxonomy that defines a uni- 
versal language, which allows agents to com- 
municate particular resource needs as well as 
offers. First simulation and optimization results 
were produced in the context of pharmacies 
(further agent models are to be included [17]) 
and urban crisis management, see Fig. 4 — as 
ABO-F is considered as an engine for smart 
crisis management it is improved and ex- 
tended continuously, where a particular focus 


? https://eu-neris.net/library/sra.html 

? https://www.bmu.de/en/topics/nuclear-safety-radiologi- 
cal-protection/nuclear-safety/response-to- 
fukushima/overview/ 
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lies on optimization and artificial intelligence. It 
is envisaged that ABO-F will be integrated into 
JRODOS. 


Figure 4 Pharmacies lacking power and medicine have 
blue icons, pharmacies lacking only medicine have yellow 
icons. All pharmacies offering medicine or an emergency 
power unit (EPU) have a red or blue dot resp. The edges 
describe a possible solution of allocating medicine and 
EPU. 


Nuclear Preparedness and Response 


Nuclear preparedness and response remain 
an important research topic due to research 
needs identified at the international level within 
the NERIS? (European Platform on prepared- 
ness for nuclear and radiological emergency 
response and recovery) SRA (Strategic Re- 
search Agenda) and at the national level with 
the many Nuclear Power Plants surrounding 
Germany? To support decision making KIT 
has coordinated and is coordinating all re- 
search activities related to the decision support 
system JRODOS^ (JAVA based Real-time De- 
cision Support). The system is able to support 
decisions about the introduction of a wide 
range of potentially useful countermeasures 


4 https://resy5.iket.kit.edu/JRODOS/ 
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(e.g., sheltering and evacuation of people, dis- 
tribution of iodine tablets, food restrictions, ag- 
ricultural countermeasures, relocation, decon- 
tamination, restoration, etc.) mitigating the 
consequences of an accident with respect to 
health, the environment, and the economy. It 
can be applied to accidental releases into the 
atmosphere and into various aquatic environ- 
ments. Appropriate interfaces exist with local 
and national radiological monitoring data, me- 
teorological measurements and forecasts, and 
for adaptation to local, regional and national 
conditions worldwide. 


Nuclear emergencies, as demonstrated with 
the Chernobyl and Fukushima disasters, do 
not stop at national borders. However, so far, 
no harmonized response was established in 
Europe or worldwide. To overcome this con- 
straint, The European Commission supports 
the installation of JRODOS in national emer- 
gency centres. Having one system operational, 
at least results are comparable, even if na- 
tional regulations may differ. 


m RODOS installation 
m RODOS installation — started 


» RODOS local users 


KIT has installed JRODOS in many countries, 
e.g. starting in 2020 in 6 West-Balkan coun- 
tries, 6 countries of the Gulf Cooperation 
Council, Armenia and lran. 


Ongoing research to improve JRODOS is re- 
lated to source term reconstruction [29] and 
ensemble modeling to capture uncertainties in 
the early phase of an emergency [27, 28]. 


Ensemble Modelling and Uncertainty 


The potential source term as well as meteoro- 
logical forecasts are the main uncertainties in 
the early phase when the release to the atmos- 
phere is expected or ongoing [26]. On the other 
hand, evacuation of the population is most ef- 
fective before the release has started. To in- 
vestigate these uncertainties, research was 
conducted within CONFIDENCE? and JRO- 
DOS has been expanded with ensemble mod- 
elling functionalities [28]. One of the questions 


2020 Installation in West Balkan countries, 
GCC countries, Armenia and Iran 


Figure 5 JRODOS installations worldwide. 


5 https://portal.iket.kit.edu/CONFIDENCE/index.php 
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that was raised but not answered was “do en- 
sembles capture the full uncertainty of a mete- 
orological forecast?”. 


To look into this question, KIT collected en- 
sembles from three different weather providers 


e  NOMADSS with 20 ensembles 

e Canadian” weather service with 21 ensem- 
bles 

e German? weather service with 41 ensem- 
bles 


at distinct times in March 2020. Sites all over 
the world were selected and ensemble calcu- 
lations have been performed. To investigate 
only the meteorological uncertainties, the 
source term was fixed. The results are pre- 
sented as percentiles, exceeding a particular 
dose value. Dark colors mean that many en- 
semble members fit therein, whereas lighter 
colors show areas where only some ensemble 
members predict exceedance of the threshold. 
The following figures show example results for 
a more constant meteorological situation and 
for one with changing wind directions. 


Nomads provider 


ha’ 


te = zd... 
Deme gu ES a 


DWD provider 


Even if preliminary, these two examples indi- 
cate that using one provider does not capture 
the full uncertainty that exist when considering 
other weather providers. Fig. 7 shows very 
similar results from all three weather providers; 
however, the meteorological conditions were 
quite stable. Fig. 6 on the other hand demon- 
strates that all three providers differ considera- 
bly, at least for the lower percentiles. Neverthe- 
less, Fig. 6 also demonstrates the added value 
of an ensemble assessment compared to a 
single forecast. As ensemble modelling is very 
time consuming — each ensemble member run 
is one project in JRODOS - a further research 
question for the following years is the number 
of ensemble members that is sufficient for op- 
erational use, considering the time constraints 
but reflecting still the existing uncertainty band. 


Canada provider 
Mk m 


Figure 6 Release in China, exceedance of iodine prophylaxis for children (B0 mSv). 


6 https://nomads.ncep.noaa.gov/ 
7 https://weather.gc.ca/canada_e.html 
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8 https://www.dwd.de/EN/Home/home node.html 
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Nomads provider 


E 


Figure 7 Release in Germany, exceedance of iodine prophylaxis for children (50 mSv). 
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Introduction 


The Hydrogen Group continued to develop 
models and in-house specialized reactive CFD 
codes (GASFLOW and COM3D) and con- 
ducted a huge experimental program to 
deepen the understanding of the mixing be- 
havior and the transient combustion phenom- 
ena, flame acceleration and deflagration-deto- 
nation-transition, of accidentally released 
hydrogen. These activities were focussed on 
phenomena related to cryogenic hydrogen in 
2020. 


Members of the group are actively transferring 
their insights and expertise into the respective 
IEA Hydrogen Task 37, the standards commit- 
tees of ISO (TC 197), CEN/CENELEC (TC 6) 
and in the European Hydrogen Safety Panel. 
Among other third party funded projects the 
Hydrogen Group has been conducting ambi- 
tious experimental programs for the European 
Fuel Cell and Hydrogen Joint Undertaking 
(FCH JU) funded projects PRESLHY - Pre-nor- 
mative research for the safe use of liquid hy- 
drogen - and HyTunnel-CS - Pre-normative re- 
search for safety of hydrogen driven vehicles 
and transport through tunnels and similar con- 
fined spaces. 


This report present some highlights from the 
experimental campaigns with liquid, cryogenic 
respectively, hydrogen conducted in the frame 
of the European PRESLHY project. The exper- 
iments will be presented with increasing con- 
finement of the reactive mixtures and start with 
a description of pool experiments. 


Liqud Hydrogen Pool Experiments 


The experimental setup for the Pro Science 
(PS)/KIT POOL experiments is shown in Fig- 
ure 1 below. The facility comprises mainly of 
an insulated square stainless-steel box (50 cm 
x 50 cm, height: 20 cm) on a scale. Three iden- 
tical boxes were made and filled up to half the 
height (10 cm) with different substrates (con- 
crete, sand, water and gravel) prior to the ex- 
periments. The concrete pool was prepared al- 
ready more than one month prior to the 
experiments, while the experiments with water 
and gravel were performed with the same box, 
which was freshly filled with the respective 
substrate for the corresponding experiment. 


LH -Trailer 


ALGU 0612017 
FR 4377 


Figure 1: POOL -facility set-up for the unignited and ig- 
nited LH2 pool experiments of the PRESLHY project. 


Inside and above the substrate, as well as 
above the facility in total 24 thermocouples 
were distributed. Above the facility also three 
measurement positions for remote H2-concen- 
tration measurements were placed, see Fig. 2. 
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Further instrumentation included optical obser- 
vation and ambient wind measurements. In 
four of the ten experiments also a large fan was 
installed to investigate the influence of side 
wind on the formation and evaporation of LH2- 
pools. 


1000 - - e 
550 om o e e 
en 
en 
250 o 
e 
100 o Pool 


Box 
(stainless steel) 
Insulation 
(PU-Foam) 


Palette 
(for transport) 
Balance 


Figure 2: Sketch of POOL-facility set-up 


In most of the experiments the pool was filled 
three times with LH2 until it started to overflow 
to investigate the influence of different initial 
temperatures of the substrate. After the fillings 
the pool was left to evaporate. The filling level 
of the pool was determined using the weight of 
the facility and the thermocouples inside the 
pool in different heights, which clearly indi- 
cated their coverage with LH2 by a constant 
value of approx. 20 K, while they started to 
show higher transient values as soon as they 
were exposed to the atmosphere, see Fig. 3. 
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Figure 3: Temperature recordings for experiment Con- 
crete02. 


With these measurements the evaporation 
rates for the different substrates have been de- 
termined, see Figure 4. Also the influence of 
side wind on this phenomenon was assessed. 
In all cases much faster evaporation rates 
were determined for side wind conditions 
(vWind & 5 m/s). A further important outcome 
of the experiments was information on the con- 
centration distribution above the pools that was 
used to define the ignition position and point in 
time for a second experimental series in which 
the gas cloud above the pools is ignited. 


Pool height from TC data 


Pool height ( 


Eats 


Figure 4: Example for determination of LH2-evaporation 
rate; here experiment Concrete02. 


With the same set-up and experimental proce- 
dure also experiments on the ignition of LH2- 
spills were performed. For these experiments 
the instrumentation above the pool was re- 
moved (compare sketch in Figure 2) and re- 
placed by 6 fast pressure sensors that were 
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positioned in heavy adapters on the ground in 
the close vicinity of the facility. The gas cloud 
above the pool was ignited using two elec- 
trodes in between which a high frequency 
spark (60 kV, 200 Hz) was generated. In total 
14 ignited POOL-experiments with the four 
substrates concrete, sand, water and gravel 
were performed. In the tests the same pools as 
in the unignited experiments were used and 
mainly ignition point in time and height were 
varied. The ignition settings were chosen on 
the basis of the results of the LH2 filling level 
measurements and the  H2-concentration 
measurements of the unignited tests per- 
formed before. 


The ignition was initiated during the evapora- 
tion phase after the second filling of the pool, 
since after the first filling vigorous boiling be- 
haviour above the still rather warm substrate 
occurs. After the second filling the boiling be- 
haviour above the precooled substrate is 
calmer and thus the filling level and concentra- 
tion measurements give more precise values. 


In the ignited experiments different degrees of 
damage were observed for the different sub- 
strate materials. For the substrates with a ra- 
ther low porosity almost no or only minor dam- 
age was observed (see left part of Figure 5), 
while the experiment Gravel04 with highly po- 
rous substrate of gravel showed a complete 
destruction of the facility (see right image of 
Figure 5). 


The reason for this exceptional behaviour, 
shown with a sequence of frames from the high 
speed video in Figure 6, is potentially con- 
densed or frozen oxygen located in the free 
space between the stones of the substrate 
layer. This mixed with the refilled LH2 and fi- 
nally reacted in the observed explosion. 


Gravel04 


Figure 5: Different degrees of damage to the facility 
observed in the ignited pool experiments for the differ- 
ent substrates 
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A i 
Figure 6: High speed video sequence of the experiment 


with ignition Gravel04 (PHOTRON High-Speed Camera, 
2000 fps) 


In summary, the observed evaporation rates 
provided valuable validation data for existing 
correlations and the observed massive explo- 
sion of the gravel substrate will initiate some 
modified recommendations for suitable ground 
materials for LH2 installations. 


Transient cryogenic jet ignition 


Transient cryogenic hydrogen jet fire behav- 
iour, including scaling and radiation properties 
has been investigated with the small scale dis- 
charge facility DISCHA. The discharged hydro- 
gen inventory varied with the initial pressure 
and temperature from 4.4 to 138 g. Similar as 
for the unignited discharge experiments done 
in 2019 [1] the release nozzle inner diameter 
was varied from 1-4 mm. Measurements con- 
sisted of background imaging system (BOS) 
combined with a high speed camera, fast pres- 
sure sensors and thermocouples. Additionally, 
a fast thermos-vision FLIR camera allowed 
monitoring the transient temperature fields. 
With mass flow rate and hydrogen distribution 
profiles known from the precursor unignited ex- 
periments, transient ignition phenomena and 
further flame development have been investi- 
gated with respect to maximum combustion 
pressure, temperature and heat flux radiation 
for model validation and hazard distance eval- 
uation. 


Unburned cold hydrogen jet was ignited with 
different delay time after jet initiation at differ- 
ent distances from the nozzle. Then, either a 
steady flame or a strong explosion with for- 
mation of a spherical shock wave might occur 
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just after ignition (Figure 7). The over-pressure 
from 0.04 to 0.115 MPa corresponds to a visi- 
ble shock wave velocity from 390 to 480 m/s 
measured by high speed BOS imaging. 


Figure 7: A shock wave formation (left) and a stationary 
jet fire (right) established under ignition of 4-mm nozzle 
and 20 MPa pressure hydrogen release: SW -shock 
wave; CH2 —unignited hydrogen 


A sequence of frames with temperature distri- 
bution is obtained by thermo-vision FLIR-Cam- 
era (15 fps) (Figure 8). For ambient conditions 
285K the local maximum combustion tempera- 
ture changes from 1100 to 540K correspond- 
ing to maximum heat flux of 85 kW/m2. At the 
same time, the average integral heat flux of 
whole surface is about 6.5 kW/m2. At cryo- 
genic temperature of 80K, the maximum tem- 
perature changes from 1330 to 710K corre- 
sponding to maximum heat flux of 177 kW/m2 
in the center of jet fire. The average integral 
heat flux of whole surface is about 11 kW/m2. 
The reason of such difference is four times 
larger hydrogen inventory and also 2.5 times 
higher mass flow rate at cryogenic temperature 
leading to 1.3 times higher temperature, 2 
times higher heat flux of flame radiation, 1.5 
times larger flame length and 1.4 times longer 
release time. 


Figure 8: Temperature distribution within jet fire structure. 
Jet fire radiation for 2-mm nozzle and 10 MPa of initial 
pressure at T=285 K (upper sequence of frames) and 
T=80 K (lower sequence of frames). 
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Using the pressure records and the results of 
extensive optical observation with up to five 
cameras the ignition and combustion behav- 
iour of the jets could be analysed. After ignition 
the flame either burns back to nozzle or not. 
This difference is important since in case the 
flame burns back to the nozzle it will continue 
burning, while in case the flame does not burn 
back to the nozzle it will quench after the igni- 
tion source is turned off. Using sensor records 
and optical observation on flame behaviour 
tests can be divided in two groups, where the 
flame burns back to the nozzle or not. When 
this data is combined with the results of the H2- 
concentration measurements of the unignited 
tests a clear trend can be observed. For very 
low H2-concentrations in the ignition position 
no ignition occurs, while for H2-concentrations 
lower than approx. 10 Vol% the jet is ignited 
but does not burn back to the nozzle. For H2- 
concentrations higher than approx. 10 Vol% in 
the ignition position the flame burns back to the 
nozzle. 
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Figure 9. Maximum measured overpressures in exper- 
iments at ambient (top) and cryogenic temperature 
(bottom) of the discharged gases. 


Concerning the combustion overpressures the 
following observations are made. In general, 
the maximum measured overpressures in- 
crease with increasing nozzle diameter, while 
the gas temperature in the reservoir seems to 
have no significant influence on the maximum 
pressure load. At ambient gas temperature the 
measured overpressures increase with de- 
creasing ignition distance (Figure 9 top), while 
in the cryogenic tests the highest overpres- 
sures are measured for an ignition distance of 
62.5 cm (Figure 9 bottom). 


CryoTube Experiments 


More than 100 experiments were made with 
the cryogenic combustion tube CRYOTUBE 
with L=5 m and 54-mm id and three blockage 
ratios BR = 0, 30 and 60% (Figure 10). About 
half of the experiments were made with hydro- 
gen-air mixtures at cryogenic temperatures 
(approx. 80 to 130 K). The temperature was 
supported by the location of the tube under a 
layer of LN2 at T=77K (Figure 11). Subsonic, 
supersonic deflagrations and detonations were 
monitored for cryogenic hydrogen combustion 
by light sensors and pressure gauges located 
along the tube. 
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Figure 10: CRYOTUBE experiments mounted on a rack 
for the ambient warm pre-cursor tests. 


Figure 11: CRYOTUBE immersed in a bath of LN2 with 
supporting structure. 


The critical conditions for flame acceleration to 
the speed of sound or detonation were evalu- 
ated as a function of initial temperature. Partic- 
ularly at 100 K, it shows a much higher hydro- 
gen concentration of 16% H2 leading to sonic 
deflagration than that of 9.6%H2 predicted by 
advanced extrapolation of existing high tem- 
perature data before the tests. The correlation 
for the critical expansion ratio based on current 
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experiments is quite simple and useful in a very 
wide range of initial temperatures (T = 90 - 650 
K): 


c*-2200T ^" 


The run-up distance to the speed of sound or 
detonation in a smooth channel (BR=0) at cry- 
ogenic temperatures is two times shorter than 
at ambient temperature. For the first time, a 
steady-state flame propagation regime in a 
smooth channel with the speed of sound in 
combustion products was frequently observed 
in cases with suppressed detonation. 


The detonation cell sizes at cryogenic temper- 
ature T = 100K have been evaluated on the ba- 
sis of existing criteria for detonation onset in 
smoot and obstructed tubes and can be ap- 
proximated by a polynomial function of hydro- 
gen concentration [H2]: 


A[mm]= 0.0006724[H2]* — 0.1039[H2]? + 
6.0786[H2]? — 159.74[H2] + 1603.3 


It appears that the detonation cell sizes for hy- 
drogen-air mixtures at cryogenic temperature 
T=100K only insignificantly differ from that at 
ambient conditions. With the correct detona- 
tion cell size, all known DDT criteria can be 
used to assess the detonability of hydrogen — 
air mixtures in different geometries and scales 
at cryogenic temperatures. 


The maximum combustion pressure at cryo- 
genic temperatures is 2-3 times higher than 
that for ambient conditions. Theoretically, the 
maximum combustion pressure for sonic def- 
lagration, which is roughly equal to adiabatic 
combustion pressure, is proportional to the 
temperature ratio TO/T, where TO =293K and T 
is the cryogenic temperature. Similarly, the 
maximum detonation pressure can also be 
predicted. It demonstrates a very high level of 
hazard under cryogenic hydrogen combustion. 
Even a sonic deflagration can be 1.5 times 
more dangerous than the detonation of the 
same hydrogen-air mixture at ambient condi- 
tions. 


Group: Hydrogen 


A broader overview on all PRESLHY results 
may be found in [2], [3] and in particular in [4]. 
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First experimental results of a supercritical Organic Rankine 
cycle using propane as working fluid 


Luciano Gardella, Hans-Joachim Wiemer 


In geothermal power generation, hot water is 
pumped at around 100-200?C from deep lay- 
ers of rock in the earth's interior and converted 
into electrical power via a cycle process, which 
is usually fed into the public power grid. Elec- 
tricity generation for these low-temperature cy- 
cles is nowadays usually realized in low to med 
enthalpy wells via an ORC process (Organic 
Rankine Cycle). 


Turbine 


Condenser 


Heat Excharger 


Fig. 1: ORC process 


The variation of the boundary conditions and 
the working fluids used has a large influence 
on the net electricity yield and thus on the eco- 
nomic efficiency of a geothermal power plant. 


The approach followed by the ITES-EVT group 
is to maximize electricity production by select- 
ing a heat transfer fluid suitable for the site- 
specific boundary conditions in Europe and a 
supercritical live steam pressure. In previous 
studies [1],[3] Vetter et al. have shown that sig- 
nificant performance improvements are possi- 
ble at supercritical conditions. 


The working fluid selection has taken into ac- 
count environmental and operational criteria, 
such as high thermal conductibility, low spe- 


cific volume, high chemical stability, low corro- 
siveness, low flammability, toxicity, low Ozone 
Depletion Potential (ODP) and Global Warm- 
ing Potential (GWP) [2]. 


In order to optimize and validate this simulation 
model we built a supercritical ORC (MoNiKa). 


The MoNiKa power plant 


MoNiKa (Modular low-temperature cycle Karls- 
ruhe) is a facility built at KIT campus north with 
the idea of studying and optimizing the ORC 
process. This installation is a small and com- 
pact power plant. It was designed as a modular 
installation to allow the study and investigation 
of different components and operational pa- 
rameter for the research of geothermal power 
generation from low-temperature heat 
sources. 


Fig. 2: MoNiKa -ORC power plant 


Components description 


The pumping system is compound by two 
pumps. The main pump is a LEWA triplex 
M514US G3G. It is a piston pump of max. 75 
kW with a maximum mass flow of 3.6 kg/s and 


59 


Group: Energy and Process Engineering 


a design pressure of 6.5 MPa (in this will be 
referenced as main pump). 


The heat exchanger manufactured by 
Gesmex, is the connection between both cy- 
cles, (thermal water and organic). It is a coun- 
ter-cross-flow heat exchanger designed to 
work in a supercritical regime. This component 
consists of 200 circular welded stainless steel 
plates grouped in five stage. The design ther- 
mal power installed is 1000 kW for full load op- 
eration 


The condenser manufactured by KUHL- 
TURMKARLSRUHE has a support sub frame 
with 3.5 m of height. It is located between the 
exit of the turbine (or throttling valve), and the 
propane tank. However, the condenser is pre- 
pared to work with water spray. The heat ex- 
change areas are built symmetrically in "V" 
configuration. They include three chambers. 
Each one is equipped with a vertical fan (im- 
peller and diffusor) of 2.5m diameter. The 
power consumption of each motor is max.13 
kW at 322 RPM, and a maximum volume flow 
rate of 44 m?/s each and in total 39 kW - 132 
m?/s. 


The last component of the circuit is the turbine 
from M+M. The axial turbine has four stages 
with a nominal rotational speed of 9996 1/min. 
The 150 MVA synchronous generator is con- 
nected to the turbine with a gearbox. The de- 
sign requirements assumed an isentropic effi- 
ciency at full load of 0.8 and a vapour quality 
limit at the turbine outlet of 0.9. 


Fig. 3:Propane tank und Condenser 
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Fig. 5: Turbine with gearbox and generator 


In addition to these main components, the fol- 
lowing auxiliary units are also important for the 
power plant operation. 
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Support pump from Grundfos, which ensures 
the upstream pressure to the main pump. 


Propane separator from Borsig for separating 
the nitrogen propane mixture from the gas 
charged in the turbine shaft seal. The sepa- 
rated propane is returned to the propane tank. 


Propane tank with a volume of 2.4 m?, also an 
additional storage tank for refilling propane. 


Furthermore, we have a 7-m fuel oil tank to 
supply the heating station. 


The control and regulation technology con- 
tainer with a redundant S7-416 PLC and the 
SIEMENS power plant software T3000. 


The Control room container with three work- 
stations for power plant and instrumentation 
monitoring. 


Scientific measuring points for the acquisition 
of additional parameters such as ambient tem- 
perature and humidity, pressure difference 
condenser, temperature and velocity field con- 
denser as well as an optical access to the pro- 
pane circuit for flow characterization. 


Instrumentation 


The sensor technology installed in MoNiKa ful- 
fills two requirements: the first for power plant 
control and regulation, the second as a plat- 
form for scientific studies conducted in the 
plant. Therefore, the plant is much more exten- 
sively instrumented than in a normal power 
plant. For balancing purposes, sensors for 
temperature and pressure are installed along 
the entire cycle, i.e., at the inlet and outlet of 
each component. 


The mass flow rate and density of the working 
fluid are measured between the outlet of the 
main pump and the inlet of the heat exchanger 
using an E+H Promass 83F sensor. This sen- 
sor uses Coriolis forces and resonant fre- 
quency to provide a direct measurement of 
mass flow, velocity and density. 


AWIKA TR34 Class A-PT100 is used to meas- 
ure temperatures at all points in the system. 
These sensors are connected as 4-wire re- 
sistance thermometers and have an operating 
range of -50 to 250 ?C. These thermometers 
are characterized by a compact design, are re- 
sistant to vibration and have a fast time re- 
sponse for operating measuring points. Type K 
thermocouples are used in areas with in- 
creased time requirements. 


For pressure measurements at the MoNiKa 
plant we use the Vegabar 81 and 82. These 
pressure transmitters are universally applica- 
ble for the measurement of gases, vapors and 
liquids. They have a ceramic measuring cell, 
which allows the sensor to behave well in cor- 
rosive and hot environments. The main differ- 
ence between purchased sensors is the tem- 
perature range in which they can operate. In 
addition, the pressure range in which they op- 
erate is different within the ORC cycle, result- 
ing in pressure measurements in the low-pres- 
sure section having higher absolute accuracy 
than in the high-pressure section. 


The sensors are in direct contact with the liquid 
(propane and water). There is no sleeve or cap 
between them. This direct contact allows better 
measurement of the fluid properties, but 
makes it impossible to remove the sensor from 
the pipeline for maintenance without first drain- 
ing it. For this purpose, MoNiKa's measure- 
ment system is designed with redundancy, 
which means that two or three sensors of the 
same type are installed at each measurement 
point. This redundant configuration makes it 
possible to increase the availability of the 
measurement technology, which on the one 
hand increases safety during operation and on 
the other hand helps to detect systematic er- 
rors. 


To measure the environmental conditions, the 
system is equipped with an EE33 series hu- 
midity and temperature transmitter. In addition, 
a set of type k thermocouples is installed in 
each chamber of the condenser. 
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Tab. 1: List of Sensors 


Magnitude Model Measure ond 
ment Range error 
Pressure — Vegabar81 | 0to 10 MPa 0,02 MPa 
Pressure Vegabar82  0to 10 MPa 0,01 Mpa 
Pressure  Vegabar82 -1to 10 MPa 0,01 Mpa 
0,005 
Pressure — Vegabar81 -1 to 2,5 MPa MPa 
0,003 
Pressure — Vegabar82 -1 to 2,5 MPa Mpa 
Promass 0,005 
mass flow 83F - kg/s 
; Promass 
density 83F 0,01 kg/L 
» TR34 
Tempera -50t0 250°C — 0,4*C 
ture class A 
m " 
[1a] m In 
a == za Ba] 
[1b] [2] 


Fig. 6: Operational measurement points at MoNiKa 


The accuracy of the data acquisition system 
is mainly affected by the analog to digital con- 
verter with a resolution of 14 bits, which corre- 
sponds to an electrical current value of +1.221 
uA. 


MoNiKa's PROFIBUS has a cycle time of 
100ms. The usual average data sampling time 
was 1 sample/ Sec. For long test runs we use 
1 sample / 10 sec. because the smaller value 
exceed the maximal file size of the standard 
data export function in T3000. 


62 


Operational design parameter 


The plant is a binary cycle, where the geother- 
mal heat source is simulated. A hot water boiler 
with one MW of thermal power heats the water 
at the site, simulating the geothermal water. 
We are able to modify the temperature and 
mass flow rate at the water loop in order to 
have a range of input conditions and provide 
different scenarios. The following table shows 
the expected design parameter at full load op- 
eration. 


Tab.2: Designed operation values. 


ORC Cycle 


Propane as Organic working fluid 


5.5 MPa 

and 
Live steam point [3] 117 °C. 
ORC mass flow 2.9 kg/s 
Turbine isentropic efficiency 80% 


Fluid quality (at the condenser in- 
let) 95% 


Pumps isentropic efficiency 70% 
Thermal Water Cycle 

mass flow 2.4 kg/s 
Temperature in 150 °C 
Temperature out 47 °C 
Pressure in 0.9MPa 
Power 

Thermal power: ~ 1000 kw 
Heat released to the ambient ~ 930 kW 
Gross Power generation ~ 150 kw 
Net Power generation ~110 kw 
Thermal efficiency ~15% 


Simulation Software GESI 


GESI (Geothermal Simulation) Vers.2.3.6b [5] 
is an in-house program that has been devel- 
oped in MATLAB by the ITES (Institute for 
Thermal Energy Technology and Safety) for 
studying and optimizing the thermodynamic 
process of the Organic Rankine Cycle (ORC). 
The software use the physical property data 
taken from REFPROP 9.0, from the National 
Institute of Standards and Technology [6]. 


This tool solve the energy and mass balance 
of the ORC Power Plant in stationary regime. 
As an input the definition of the thermal water 
values, selection of the operational points, the 
ambient characteristic and the equipment is 
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necessary. The tool provides as result the ma- 
jor values of the whole process and deliver di- 
agrams e.g. temperature- and entropy dia- 
gram. 


The results of the simulation allows to validate 
the power plant model and to determine the ef- 
ficiency ofthe pump and turbine in full- and part 
load operation. 


Test run 


MoNiKa is a research platform that was com- 
pleted at the end of 2018. The first tests were 
performed in November 2019, without the tur- 
bine, operating in bypass. The focus of these 
tests was to investigate the entire cycle and all 
the components involved. 


Despite the fact that the plant was fully com- 
pleted, all components are installed and func- 
tional, some control loops of the control system 
were not fully optimized and automated at this 
time. This means that most of the test runs 
have been manually controlled. 


With the turbine tests in January 2021, the en- 
tire commissioning of the plant was completed. 


The optimization of the control and regulation 
circuits is a task in parallel to the actual test 
operation of the plant and will be further devel- 
oped on an ongoing basis. 


The first longer bypass test took place on 
23.01.2020 and lasted 20 hours with an ORC 
mass flow of 10096,7096 and 50 %. The table 3 
shows the results for the 10096 case. 


Table 3: test run results 10096 load 


Full 
Load Accuracy Unit 
100%m 
ORC 
cycle 
m onc 2.87 +0.03 kg/s 
T AIR cond 
out 6.79 10.75 °C 
T[1] 13.8 10.2 °C 
p[1] 0.722 10.003 MPa 
T[1a] 14.3 10.2 °C 
p[1a] 0.721 10.003 MPa 
T[1b] 14.0 10.2 °C 
p[1b] 0.777 +0.003 MPa 
T[2] 17.0 10.2 °C 
p[2] 5.53 +0.01 MPa 
T[3a] 108.1 10.4 °C 
p[3a] 5.52 10.02 MPa 
T[3b] 108.3 10.4 °C 
p[3b] 5.50 +0.02 MPa 
T[4] 14.6 10.2 °C 
p[4] 0.72 +0.01 MPa 
Water 
Cycle 
m w 2.4 +0.02 kg/s 
Piw 0.9 +0.01 kg/L 
T win 150.2 +0.5 °C 
T wout 61.74 +0.3 
Ambient 
Tamb 0.98 +0.75 °C 
Main 
pump 
Phya 26 +2 kw 
Peiec 27.5 +0.2 kW 
Support 
pump 
Phya 0.31 +0.006 kW 
Peiec 1 .06 +0.02 kW 
Heat 
transfer 
He Ex 894 +8 kW 


With these measured data we can calculate 
the enthalpy and entropy values for the four cy- 
cle points as mentioned in Fig.1. 


63 


Group: Energy and Process Engineering 


Table 4: Test run enthalpy and entropy results with the 
estimation for point 4. 


ORC 


points Enthalpy Entropy Quality 
kJ/kg kJ/kg K - 
[1] 235 1.12 Subcooled 
[2] 245 1.13 Subcooled 
Super- 
[3] 557 2.04 heated 
[4] = - - 
Aa 557 2.24 0.925 


At point 4 the bypass outlet/condenser inlet we 
have only pressure and temperature sensors. 
In the two phase region we need an additional 
parameter to define the enthalpy or entropy 
values. 


So we can estimate the enthalpy with the as- 
sumption of an isenthalpic throttling from point 
3 to 4 (bypass valve) and neglect the heat loss. 


The heat transfer released from the water cy- 
cle is determined to 899 kW. The amount of 
heat absorbed by the propane is 886 kW with 
a deviation of 10 kW. This mean that we have 
a heat transfer efficiency of 98% +1.5%. 


Nevertheless, the total heat transfer is about 
10% less than the design value. This circum- 
stance significantly affects the performance of 
the entire power plant. It is not only the power 
generation that is affected. The biggest impact 
is the quality of the fluid after expansion. 


The pressure loss of the heat exchanger is less 
than the accuracy of the pressure sensors. 


The designed minimal temperature difference 
of the heat exchanger was 10 K and we meas- 
ure a value of 33.5 K. 


The total pump efficiency (hydraulic power / 
electric power) for the main pump was above 
90% in the range of 50% -100% of ORC mass 
flow with a maximum value of 97%. 


For the support pump, we determine the total 
efficiency to 29% -33%. 
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In order to balance the released heat in the 
condenser, the heat release on the airside was 
verified. 


Using the air temperature measurement and 
under the boundary condition that the fans 
were operating at full power, the heat released 
was calculated and then the thermodynamic 
values of the state point 4. 


The result is inconsistent, since from the simu- 
lations a heat release of 1100 kW was ex- 
pected and in the test run we measured only 
about 930 kW as maximum released power. 


These results show a similar steam quality af- 
ter throttling and a slight sub cooling of the pro- 
pane at the condenser outlet, which is not con- 
sidered in the simulation. However, the 
difference in enthalpy and entropy is only 0.196 
and does not explain the lower heat release. 


Additional investigations have shown that the 
propane tank, due to the large propane mass, 
determines the condensation pressure at the 
condenser outlet and thus also the minimum 
temperature. 


More detailed investigations of the condenser 
will be carried out in the future. 
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Fig. 7: Ts-Diagram of the ORC in bypass operation in 
comparison to the designed bypass operation. 


During turbine commissioning in December 
2020, significant vibration amplitudes occurred 
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in the turbine and gearbox bearings, exceed 
the predefined limits and causing the turbine to 


trip. 


Only after readjustment and replacement of 
the vibration sensors, it was possible to run the 
turbine up to nominal speed. 


Unfortunately, the synchronization of the gen- 
erator only worked after a manual adjustment 
of the excitation current controller, so that we 
were able to feed electricity into the grid for the 
first time on January 25, 2021. 


In this test of commissioning, the turbine was 
loaded up to an electrical gross power output 
of 88.6 kW. 


Summary 


The MoNiKa cycle in stationary bypass opera- 
tion was studied. As a first step, we study the 
installed sensors in the facility in order to de- 
termine the accuracy of the measurements. 


Experiments were carried out at three operat- 
ing points 50%, 70% and 10096 of ORC mass 
flow and analyzing the behavior of the entire 
power plant as well as the individual compo- 
nents. 


The 1096 less heat transfer at the heat ex- 
changer in comparison to the design values 
makes it necessary to redefine the optimal full 
load point and limit the range of operation due 
to the minimal allowable vapor quality at the 
turbine in- and outlet. 


The adaption of ORC turbine inlet live steam 
conditions with an ORC mass flow of 2.7 kg/s, 
which is only 9396 of the designed mass flow, 
is the result of the simulations performed with 
GESI with the corrections already discussed, 
and considering the limits of the thermal water 
cycle. 


We completed the commissioning of the 
Monika plant with the successful power grid 
supply. In the further experimental operation, 


the focus will be on investigations and optimi- 
zation of the single components and the full 
system performance. 
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Annual report of the Framatome Professional School 


Andreas Class 


The Framatome Professional School (FPS) or- 
ganizes training courses in safety technology, 
thermal hydraulics, reactor physics, stress 
analysis, and hydrogen technology etc. 
Courses are tailored for engineers specializing 
in new topics within Framatome and they are 
open for students from KIT and other universi- 
ties. Moreover, through the regular exchange 
between KIT and Framatome public funding 
opportunities for joint proposals are identified. 


The FPS coordinates the InnoEnergy Master 
Schools Programmme "Energy Technologies" 
which offers a broad-based education in a va- 
riety of key engineering disciplines. A focus 
lays on the inter-disciplinary aspects of the 
sustainable energy sector and an in-depth un- 
derstanding of the role of innovation and entre- 
preneurship in the future energy industry. Stu- 
dents receive a Master's degree from two of 
the four participating universities (IST: Instituto 
Superior Técnico in Portugal, Grenoble INP: 
Institute of Technology in France, UU: Uppsala 
University in Sweden, or KIT). The FPS also 
organizes ex-curricular activities for students 
such as the annual SiemensEnergy New En- 
ergy Challenge. 


The FPS has a long history of applied research 
and strong links to other ITES-groups exist. 
The FPS applies Computational Fluid Dynam- 
ics and System Codes, develops Reduced Or- 
der Models for Uncertainty Quantification and 
mesh-free modelling of three-dimensional ther- 
mohydraulic components, and exploits meth- 
ods of applied mathematics such as asymp- 
totic series expansions to complex physical 
problems. 


The applied focus of the research of the FPS 
led to inventions in 2020. A lab-on-a-chip 
STOKES II developed by P. Marthaler and A. 


Class allows driving um-scale fluid flow on a 
chip in a freely programmable fashion. This 
provides a new platform technology for com- 
plex diagnostics and chemical synthesis tasks. 


The FPS invention BattMarines addresses 
electric energy storage in a power-heat-power 
concept. 


In the following sections a number of activities 
of the FPS are summarized. 


CFD Validation of heated rod bundles 
in heavy liquid metals 


Abdalla Batta, Karsten Litfin, Andreas Class 


In our group, validation of CFD models for ther- 
mal-hydraulic simulation is an essential part of 
our work. Recently many studies are per- 
formed for supporting the design of next gen- 
eration reactors like the "Multi-purpose hybrid 
research reactor for high-tech applications 
(MYRRHA) planned at SCK-CEN, Belgium. 
For this purpose a series of EU-projects have 
been started in the last years like EURO- 
TRANS (2005-2010), THINS (2010-2015), 
SEARCH (2011-2015) and MAXSIMA (2012- 
2018), to name a few. 


The ongoing Project PATRICIA (2020-2024) 
aims to support safety studies for heavy liquid 
metal cooled fast reactor systems. The effect 
of a well-defined porous blockage in a wire 
spaced 19-pin rod bundle and the effect of rod 
bundle deformation will be experimentally in- 
vestigated at KALLA with numerical support 
from our group. Further numerical support will 
given for the project PASCAL (2020-2024) 
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Figure 1: Schematic side view of a rod bundle test section in heavy liquid metal with two small blockage elements. 
All axial positions are expressed in mm, referring to the onset of the heated zone (z=0). 


where the effect of deformation on the local 
flow and temperature field will be studied. 


An example for the validation of experimental 
results obtained within the EU FP7 project 
"Methodology, Analysis and eXperiments for 
the Safety In MYRRHA Assessment 
(MAXSIMA) project will be presented in the fol- 
lowing subsection. More details about the pro- 
ject can be found elsewhere [1], [2], [3]. Figure 
1 shows the test section with two blockage el- 
ements (marked in yellow) installed at selected 
locations in the heated section (marked in red). 
Details of the study are given in [4]. 


Model validation 


In this example, we present a comparison of 
the experimental results to CFD calculations. 
To quantify the impact of the blockage (I) onto 
the local temperature we define the impact fac- 
tor as: 
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Tpiockea = Tunblocked 
I = 100 —— 
Tbtocked — Linlet 


The impact factor is presented for both the ex- 
perimental and numerical data in Table 1. The 
evaluation employs temperature values meas- 
ured behind the first blockage (Central one). 
From experimental data, one can see that wall 
temperature is influenced with maximum 
10.9 % at the thermocouple Twi2 R16 24o, where 
the blockage impact on fluid reaches 22.7 % at 
TmL2_sa1. Considering the numerical results, it 
can be seen that the blockage impact is 
smaller than the experimental values and take 
place at different location. This shows that cur- 
rent models are not yet adequate to predict 
flow mixing around blockages. 


For the comparison between numerical and 
experimental local results, the error (Err) can 
be calculated by: 


Texp erimental T Tnumerical 


Err = 100 


Texp erimental ^ *inlet 
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Table 1: Experimental and numerical temperature data for case 157 and 24 of the blocked rod bundle experiment at 


KALLA. All temperatures in [°C] 


Experimental Numerical Num. error comp. to Exp. 
type S TR l Pe] Hine T dd 3 b 
Case number 157 24 157 24 157 24 
| ML2 R1 270 278.31 278.64 -0.4 285.7 283.6 2.5 9.5 6.3 
| ML2 R6 330 276.64 276.56 0.1 281.6 281.0 0.7 6.4 5.8 
| ML2 R7 210 269.04 269.67 -0.9 274.7 278.7 -5.0 8.2 12.9 
| ML2 R16 240 257.11 251.51 10.9 246.9 243.1 8.7 -17.9 -16.2 
| ML2 R16 330 266.81 267.83 -1.5 264.1 268.5 -6.4 -4.1 0.9 
| ML2 R17 30 264.27 266.41 -3.2 269.4 271.7 -3.3 8.0 8.0 
| ML2 R18 150 263.95 266.2 -3.4 276.4 274.3 2.8 19.4 12.2 
| ML2 S5 262.83 260.75 3.4 274.4 266.3 12.3 18.5 9.1 
| ML2 S20 270.01 259.82 17.0 262.0 262.7 -1.1 -11.4 4.8 
|I ML2 S21 261.59 257.19 7.7 262.0 258.9 5.2 0.7 3.0 
| ML2 S33 245.64 246.52 -1.9 259.3 249.9 18.9 30.0 7.3 
[T_ML2_S41 237.57 230.61 22.7 225.8 223.3 10.7 -31.2 -23.8 


The calculated error values are given in the last 
two columns of the Table. Note that the maxi- 
mum error is observed for the fluid tempera- 
tures measured near the outer wall. This can - 
to some extent — be related to heat losses at 
the outer wall of the bundle which has been ne- 
glected in the simulation. However, the domi- 
nant reason is a weak capability of the model 
to predict turbulent liquid-metal flow with obsta- 
cles and strong circulation. 


Based on the presented data the Nusselt num- 
ber can be calculate yielding for the experi- 
mental data: Nubtockea = 10.54 and 
Nuunblockea = 11.03. The value for Nu,nptocked 
is similar to the value Nu = 11.2 calculated 
based on comprehensive data. From the nu- 
merical data Nuptocked = 9.73 and 
Nuunblockea = 9.88 are obtained. Comparing 
experimental to numerical results for Nusselt 
numbers better agreement than for local tem- 
peratures is found. Thus, the blockage has a 
minor effect on the calculated Nusselt number, 
since the latter is based on the average tem- 
perature. Note that the predicted temperature 
error scatters exhibiting negative and positive 
values. Therefore, in the averaged values the 


difference relative to simulated results is less 
when comparing local data to experimental 
data. Detailed discussion on the experimental 
results can be found in [5] 
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Creation of a reduced order model for 
the k - model 


Jorge Yanez, Andreas Class 


Introduction 


Nuclear technology is a sector in which safety 
is insured in a very strict manner. In these 
frames, is becomes crucial to analyze the sta- 
tistical trustworthiness of Computational Fluid 
Dynamics (CFD) results. That is to assess the 
effect of the bias and the unavoidable random 
variation and aleatory disposition of the initial 
conditions of a realistic scenario. For this task 
Reduced Order Models (ROM) come to hand 
[3]. Here we address the creation of a ROM for 
the standard k-e- turbulence model as it repre- 
sents the “work horse” of CFD analysis. This is 
not a simple task since the k-e model equa- 
tions exhibit multiple non-linearities. 


Theoretical model 


We start with the formulation of the Navier- 
Stokes equation in its Boussinesq approxima- 
tion and the complementary k-e equations[4]. 
In our methodology [2], we utilize the Method 
of Snapshots [1] to create a simplified model. 
That is, we obtain the transient solution of our 
equations with the high fidelity solver STAR- 
COM at times {t,,....,t„} . We form with them a 
snapshot matrix Y. We carry out its Singular 
Value Decomposition, Y = VZGT. We project 
the n dimensional vectorial space in a-re- 
duced-sub-space of dimension n — j dropping 
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the last j vectors in V. The first n — j vectors of 


V create a vectorial space at most Sr, Il o; 11? 


distant form the original one. 


To develop the reduced order model, we write 
the k-e equations in terms of the reduced order 
basis. We follow the same rationale as [3]. Af- 
ter a quite complex algebra, the model adopts 
the form, 


— AtG — 0 


Note that D, E, F, G and J are matrixes with a 
much smaller dimension compared to the de- 
grees of freedom of the high fidelity calculation. 
They depend on the dimension of the reduced 
basis. Yet, the matrix entries depend on the 
comprehensive high-fidelity calculation. All 
matrix entries are calculated once during the 
post processing of the high-fidelity results, and 
subsequently stored. 


Application 


We apply our construct to a convection case of 
domain 10x10x1 meter and resolution 0.1 me- 
ter. We consider the volume completely iso- 
lated for the whole set of variables, with zero 
gradient perpendicular to the boundaries. 


Considering the snapshots generated and a 
maximum distance of 1077 to the projection 
space, the dimensions of the reduced space is 
25 considering all variables of the system. 
Note that that means that the--full--matrices of 
the model are in fact of size 25 x 25 in total, 
compared with the dimension 
100.000x100.000 of the—sparse--matrices of 
the original calculation for each variable. 


The model has been implemented in an in 
house Python code. The results of the integra- 
tion of the ROM compared with the original 
CFD results can be seen in Figure 2. The CFD 
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result has been projected into the reduced ba- 
sis. The amplitudes in this basis of the ROM 
and the CFD can be directly compared. We 
see an excellent agreement between CFD and 
ROM. 
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Figure 2 Comparing the reduced order model (ROM) 
with original CFD calculations show an excellent 
agreement. (top: velocity; bottom: k) 


Conclusions 


We have shortly summarized the procedure 
we have followed for the creation of a reduced 
order model of the k-e turbulence model based 
on Proper Orthogonal Decomposition and Ga- 
lerkin projection. The model finally adopts a 
simple expression. 


In the whole derivation procedure, on-line op- 
erations, that is, those to be done during the 
integration of the reduced model, have the 
complexity of the reduced basis. Only off-line 
operations, done once during the post pro- 
cessing the results of CFD have the complexity 
of the original basis. 


The construct obtained shows a very good per- 
formance reproducing the calculation pro- 
posed to illustrate its capabilities. 
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A free programmable lab-on-a-chip 


Philipp Marthaler, Andreas Class 


Background 


Lab-on-a-chip systems are utilized in many 
fields for rapid and reproducible chemical anal- 
yses on a microscale. Especially in medical 
technology, such tools support diagnoses and 
research. The design of most systems, how- 
ever, is customized for one specific applica- 
tion. Important application cases are point-of- 
care diagnostics or DNA sequencing. Such 
single-application systems are characterized 
by networks of micro-channels. Channel net- 
works ranging from rudimentary (single-use 
systems) to complex (covering multiple use 
cases) are available. 


In the described systems, probes are embed- 
ded and transported in a fluid environment. 
The carrier fluid which contains the diffuse 
samples is excited to flow by induced-charge 
electro-osmosis. Most common is the excita- 
tion with AC potentials applied via electrodes 
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inside the channel wall. In contrast to the appli- 
cation of DC fields, that method has several 
advantages. In particular, it is applicable in 
more complex geometries and creates more 
continuous flow patterns. 


Innovative approach 


We developed a novel lab-on-a-chip system 
that does not depend on channel networks. 
The carrier fluid with diffuse probes is trans- 
ported using AC electro-osmosis with traveling 
waves, also known as traveling wave electro- 
osmosis. 


The central element of the development of that 
system is a spectral code that simulates the 
underlying electro-hydrodynamic phenomena. 
With that code, an earlier ITES patent [1] was 
examined for its effectiveness and feasibility. 
Based on the results, the excitation method 
was improved from alternating current to trav- 
eling waves. The code was further utilized to 
find the optimal operating parameters for the 
system based on fluid properties, ion species, 
and concentration. Compared to most availa- 
ble systems, high transport velocities were 
reached. 


The optimized excitation method can now be 
applied in a redesigned 1mm x 1mm chip 
(schematically shown in Figure 3) with a free 
flooded surface. Channels are not necessary 
for this system. Thus, many processes can be 
performed on one chip. Those processes de- 
pend mainly on the software implementation 
and only to a minor part on the chip geometry. 
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Figure 3 Schematic sketch of the chip structure. Dark 
blue paths represent example flow patterns in the 
flooded surface (light blue). The electrodes are con- 
trolled by single control units (black squares) and are 
responsible for flow excitation. 


Besides flow excitation and probe transport, 
the design of a wire network inside the chip 
supplying a large number of electrodes was an 
immense challenge. For the prototype design 
and numerical verification, Philipp Marthaler 
and Andreas Class won the 2nd prize at the 
KIT innovation contest NEULAND [2]. The 
award ceremony is shown in Figure 4. 


CIT 


NEULAND 


A 


Figure 4 Award ceremony of the NEULAND contest. 
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BattMarines 


Abdalla Batta, Andrea Bellelli, Karsten Litfin, 
Andreas Class 


A mayor challenge of the Energiewende is en- 
ergy storage to bridge periods where neither 
wind nor photovoltaics can provide renewable 
energy sources. For the success of storage so- 
lutions in the market, both a low price and an 
acceptable efficiency are equally relevant fac- 
tors. Power-Heat-Power solutions for energy 
storage accept compromises with respect to 
the achievable efficiency and aim at low sys- 
tem prices to become competitive. The novelty 
of the proposed BattMarines is to use a con- 
ventional Rankine-cycle in batch operation as 
an energy storage system and to exploit scal- 
ing to large dimensions (storage capacity -1 
GWd) to reach economical goals. 


Figure 5 shows a schematic sketch of BattMa- 
rines [1]. A high-efficiency steam turbine power 
plant is installed on a floating platform on the 
sea. A huge tank resides on the seabed and 
contains large quantities of hot live steam. To 
produce electric power on demand steam is 
drawn through the vertical steam line connect- 
ing the platform and the storage vessel. Charg- 
ing the system is accomplished by electrically 
heating fresh water when a surplus of renewa- 
ble energy becomes available. Even though 


the concept appears to be naive, it offers re- 
markable advantages with respect to its scala- 
bility to the anticipated large storage capacity. 
The steam within the thin-shell storage vessel 
is stored at the ambient pressure of the sur- 
rounding water at large depth and thus exhibits 
high energy density. The heat losses of the 
storage vessel and the steam line are obvi- 
ously quiet large but in comparison to the huge 
energy contents, they become acceptable, so 
that storage times ranging from several hours 
to seasonal storage become feasible. 


Figure 5 Sketch of BattMarines Power-Heat-Power stor- 
age concept. 
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